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Abstract 

The analytic contribution to the low energy expansion of Type II string amplitudes 
at genus-one is a power series in space-time derivatives with coefficients that are de¬ 
termined by integrals of modular functions over the complex structure modulus of the 
world-sheet torus. These modular functions are associated with world-sheet vacuum 
Feynman diagrams and given by multiple sums over the discrete momenta on the torus. 
In this paper we exhibit exact differential and algebraic relations for a certain infinite 
class of such modular functions by showing that they satisfy Laplace eigenvalue equa¬ 
tions with inhomogeneous terms that are polynomial in non-holomorphic Eisenstein 
series. Furthermore, we argue that the set of modular functions that contribute to the 
coefficients of interactions up to order are linear sums of functions in this class 

and quadratic polynomials in Eisenstein series and odd Riemann zeta values. Integra¬ 
tion over the complex structure results in coefficients of the low energy expansion that 
are rational numbers multiplying monomials in odd Riemann zeta values. 
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1 Introduction 


The low energy expansion of the Type II superstring effective action generates an inhnite 
set of higher derivative interactions that generalize the Einstein-Hilbert action of classical 
general relativity. The dependence of these interactions on the vacuum expectation values of 
the scalar helds is highly constrained by space-time supersymmetry and duality. The detailed 
study of these higher derivative interactions is motivated by the desire to obtain a better 
understanding of the full non-perturbative structure of string theory, and has revealed some 
striking systematics of the modular structure with respect to the dualities of string theory. 
The simplest nontrivial example is provided by uncompactihed Type IIB theory, which has 
maximal supersymmetry and for which the S-duality group is 5'L(2,Z), and the higher 
derivative interactions depend upon the complex axion-dilaton through non-holomorphic 
modular forms. In this case the low order terms in the low energy expansion preserve a 
fraction of the supersymmetry. These BPS interactions are expected to be “protected” in 
some sense, and are known to possess a small number of perturbative contributions. However, 
little is known about the structure of the modular forms associated with higher order terms. 
Important information concerning these higher derivative interactions would be provided by 
determining their contribution order by order in string perturbation theory. 

The low energy expansion of string scattering amplitudes, from which the higher deriva¬ 
tive interactions are obtained, has a rich mathematical structure already in uncompactihed 
string theory, at any order in perturbative theory. For example, the higher derivative inter¬ 
actions that arise at tree-level in Type II string amplitudes have coefficients that are rational 
multiples of odd multi-zeta values, and these in turn are related to multiple poly-logarithms. 
Contributions from higher loop amplitudes are more difficult to analyze due, in part, to the 
presence of non-analytic behavior associated with massless intermediate states. However, the 
procedure for isolating those non-analyticities and extracting local effective higher derivative 
interactions is by now well-established, at least at low orders in the low energy expansion. 

Perturbative calculations in Type H string theory in ten-dimensional Minkowski space- 
time have been carried out to various string loop orders and various orders in and have 
been compared against the predictions of space-time supersymmetry and duality when avail¬ 
able [1-12]. Specihcally, the strength of the interaction was evaluated at one-loop order 
in [1], and the cancellation at two loop order, predicted from space-time supersymmetry 
and duality, was verihed in [7]. The strength of the tree-level and two-loop orders of the 
D^TZ'^ interaction were predicted on the basis of supersymmmetry and duality in [5,6], and 
its one-loop contribution was predicted to vanish, which was verihed in [4]. The non-zero 
predicted value at two-loop order was verihed in [8]. The predictions of supersymmetry and 
duality made for the D^TZ'^ interactions in [9], were verihed at two-loop order in [12], and at 
three-loop order in [11] (although in this last instance a suspicious factor of 3 still needs to be 
elucidated). These comparisons, made originally for the case of ten-dimensional Minkowski 
space-time, have been extended to toroidal compactihcations, which are associated with 
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higher rank duality groups, in a number of papers, including [13-23]. 

One theme that enters in many of these works is the consideration of homogeneous and 
inhomogeneous Laplace-eigenvalue equations in which the Laplace-Beltrami operator acts 
either on the scalar field vacuum expectation moduli or on the moduli of the worldsheet 
Riemann surfaces of string perturbation theory. In particular, in the matching of the 
interaction at two-loop order, the relation between the integrand of its coefficient with the 
Zhang-Kawazumi invariant [24] could be used fruitfully only once it was realized that this 
invariant satisfies a Laplace-eigenvalue equation on the moduli space of genus-two Riemann 
surfaces [12]. Therefore, there appears to be strong evidence that a rich mathematical 
structure emerges not only for the low energy expansion of string scattering amplitudes, but 
also for their corresponding integrands on the moduli space of compact Riemann surfaces. 

The study of the genus-one expansion at low orders was initiated in [4,25]. The results to 
be presented here will give a more complete and unified understanding for certain classes of 
the expansion coefficients,^ and they hint at an underlying elegant mathematical structure 
that generalizes the multiple zeta values that arise at tree level [26-31]. The expansion 
of the integrand of the genus-one superstring four-graviton amplitude naturally leads to 
generalizations of the classical Kronecker-Eisenstein series [32] closely related to series that 
have been considered in [33,34] for example. We note that there has also been some recent 
progress in understanding the mathematical structure of iV-particle genus-one open string 
amplitudes [35] based on elliptic multiple-polylogarithms constructed in [36]. 


1.1 Overview of results 

In the present paper, we shall concentrate on the case of genus-one and consider a world- 
sheet torus S with complex structure modulus r. The modular functions of interest are 
then generated by vacuum Feynman diagrams for a free massless scalar field on S, and may 
be represented by multiple sums over the discrete momenta on the torus, subject to overall 
momentum conservation. First, we shall consider the infinite class of contributions arising 
from SL{2, Z)-invariant modular functions defined by the following multiple sums. 




E 

(mr,nr)^(0fl) 


^m,0 ^n,0 

r=l 


T2 


Tl\mrT + Ur-p 


( 1 . 1 ) 


where t = ti + iT 2 with ti,T 2 real and T 2 > 0. For r = 1, 2,..., p, the integers Ur ^ 
represent the world-sheet momenta p,, = rUrT -|- Ur on the torus. Total momentum conserva¬ 
tion is enforced by the Kronecker S factors which require the vanishing of m = mi 
and n = ni + ... + Up. The functions may be defined for all complex values 

of ttr by analytic continuation, though the cases of interest here will be limited to integer 


^They will also correct a number of numerical errors in earlier papers. 
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values of ar such that the inhnite momentum sums are convergent. The function 
evaluates a (p — l)-loop Feynman diagram on a world-sheet torus of modulus r. The set of 
functions Cai,...,ap dehnes sub-classes of Feynman diagrams contributing to string amplitudes 
with arbitrary numbers of external massless states. 

The weight w = ai + ■■■ + ap of the diagram corresponds to the number of scalar Green 
functions in the diagram, and labels the order in the low energy expansion at which the 
modular function contributes, namely schematically Such a diagram can contribute 

to the low energy expansion of the iV-graviton amplitude with N > w — 1. It is impor¬ 
tant to stress that the notion of weight w introduced here is different from the notion of 
modular weight under SL{2, Z) transformations; the modular weight of every G-function 
dehned above is zero. In the present paper we shall determine all higher derivative effective 
interactions up to weight w = 5 corresponding schematically to 

There are many other contributions to the coefficients of general higher derivative inter¬ 
actions that are dehned by more general Feynman diagrams. However, we will hnd relations 
between such contributions and the question of how many diagrams of weight w are gen¬ 
uinely distinct is a fascinating one. In particular, we shall examine one weight 5 contribution 
to the four-graviton amplitude which is not of the type of the G-functions above, but does 
contribute to the effective interaction 

We shall now provide a more detailed description of the results for the modular functions 
Cai, .,ap('r) obtained in this paper for increasing values of p. 

• For p = 1, the function G vanishes identically. 

• For p = 2 we have the relation, 

Ca,,a2{r) = Ea,+a2{r) ( 1 . 2 ) 

where, 

= E .1 "'1 12 . ( 1 - 3 ) 

TT^lmT + nr 

(m,n)^(0,0) 

is the SL{2, Z)-invariant non-holomorphic Eisenstein series.^ The Eisenstein series satishes 
the following Laplace-eigenvalue equation, 

AEs = s{s - l)Es (1.4) 

where the Laplace-Beltrami operator on the upper half plane is dehned by A = Ar^drdf. 
The Eisenstein series is known to be the only non-holomorphic modular function to satisfy 
(1.3) and to have polynomial growth as T 2 —)■ cxd. The Eisenstein series Eg must be assigned 

^The normalization factor of 7r“® in (1.3) is slightly unconventional but has been included for later 
convenience. Henceforth, the dependence on t will generally not be exhibited explicitly. 
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a weight w = s in line with the dehnition of the weight of C'ai,...,ap given earlier (and again, 
w should not be confused with the modular weight, which is zero for Eg). 

• For p > 3, the functions C are special cases of Kronecker-Eisenstein series [32], The 
simplest inhnite class has p = 3 for which we shall derive the most complete results in this 
paper. In particular, we shall show that each modular function Ca,b,c is a linear combination, 
with rational coefficients, of a basis of modular functions satisfying the inhomogeneous 

Laplace-eigenvalue equation, 

(A - s(s - 1)) (1.5) 

Here, w = a + & + c is the weight; s is given by s = tc — 2m where m may take one of 
the following values m = 1, 2 ,..., [(tc — l)/ 2 ]; and the multiplicity index p takes values 
p = 0,..., [(s — l)/3]. The function (^^;s;p(Es/) is a polynomial of degree 2 and weight w in 
the Eisenstein series Eg' with the integer s' in the range 2 < s' < w. 

The cases of low weight w may be written down simply in terms of the original functions 
Ca,b,c, and we will hud, for example,^ 



= 6E3 


~ 2)C'2qq 

= 9E4 - El 


AC'2,2,1 

00 

(1.6) 


For these low eigenvalues, no degeneracy occurs and the multiplicities are 1. The equations 
for higher weights are more involved and will be spelled out in section 3.5. 

The structure of the spectrum of the Laplace-eigenvalue equations, summarized above, 
will be proven with the help of a generating function for the modular functions Ca,b,c- This 
function is closely related to the sunset diagram of a three scalar helds with three different 
mass parameters, and exhibits remarkable properties with respect to the non-compact group 
SO{2, 1; M). Its structure, which is governed by the representation theory of SO{2, 1; R) mod 
three, will provide a complete explanation of the spectrum announced in (1.5). 

The Laplace-eigenvalue equation for eigenvalue s(s — 1) = 0 at each odd weight w is 
special in that it exhibits a single inhomogeneous term proportional to E^. Therefore, it 
may be integrated with the help of (1.4) and, for each odd weight w, will lead to one algebraic 
relation between the functions Ca^^c and Eisenstein series. The Laplace equation determines 
the algebraic equation up to an additive constant of integration which may be obtained from 
the asymptotic expansions of the equation. In the cases cited above, we hud the relations, 

Cl,1,1 = E3 + (^(3) 

C 221 = -^5 + ^ (1.7) 

5 30 ^ ^ 

^The structure of these equations, and the equations for higher weight functions, is analogous to that of 
equations that arise in considerations of S-duality in [9,37], Such equations have solutions that were discussed 
in [10]. 
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The result on the first line had been obtained previously by Don Zagier [38] using an artful 
direct evaluation of the multiple sums. For all higher odd weights, an explicit form may 
be obtained for the algebraic equation which results from integrating the zero eigenvalue 
equation. It will be presented in equation (3.57) in section 3. 

When s > 2 the inhomogeneous term in (1.5) is nonlinear in E^. In such cases the inho¬ 
mogeneous term in the equation for contains quadratic expressions in the Eisenstein 

series, as in the second equation of (1.6), and the solutions are no longer Eisenstein series. 
Whether the solutions to such equations can be expressed in terms of Poincare series by 
applying the methods used in [10] remains to be established. 

• For p > 4, the situation is much less well-understood, and we shall not attempt to 
present a systematic treatment here. Instead, we shall limit attention to those modular 
functions that enter into the evaluation of higher derivative interactions resulting from the 
four-graviton amplitude up to order In this case, a general diagram has 4 vertices 

and 6 edges labelled by their contribution to the weight ii,... ,ie. We shall use the notation 
introduced in [37] in which such a diagram is denoted by the symbol with 

the understanding that we drop any entry corresponding to ir = 0. In [37] various properties 
of these D functions were studied, including their asymptotic behavior near the cusp. Many, 
but not all, of the D-functions are of the form (1.1)^. 

To order which corresponds to weight 4, the non-trivial contributions are from 

-02,1,1 = ^ 2 , 1,1 and D 4 = The hrst has already been discussed. The second will be 

conjectured to satisfy an inhomogeneous Laplace-eigenvalue equation with eigenvalue 2, 

(A - 2) (5C'i, 1 , 1,1 - 15El - I 8 E 4 ) = -120 E^ (1.8) 

Evidence for the conjecture is obtained by verifying the equation in the asymptotic expansion 
near the cusp where T 2 ^ 00 , and we have done so for all pure power behaved terms, as well 
as for the leading exponential corrections of order It would be valuable to obtain a 

full proof of the conjecture, but this is beyond the scope of this paper. Equation (1.8) by 
itself does allow for a simple integration, and therefore does not imply any simple algebraic 
relation of the type that we had derived for eigenvalue 0. However, the inhomogeneous 
term that occurs in (1.8) is precisely the one that also occurs in the middle equation in (1.6), 
and upon elimination of this term, we obtain a Laplace-eigenvalue equation of the form, 

(A — 2) (C*i, 1 , 1,1 — 24C*2,i,i — 3ii/2 -|- 18774^ = 0 (1-9) 

The kernel of (A — 2) on modular functions with polynomial growth is spanned by E 2 . 
Verihcation of the power asymptotic behavior near the cusp shows that such addition is 
absent, and we obtain, 

(Fi, 1 , 1,1 = 24 C' 2 , 1,1 + 3 - 18 E 4 (1.10) 

"‘However, contrary to the assertion in [37], in order to evaluate the integral of these functions over the 
complex structure we need to determine the full non-perturbative structure of the modular functions. 



In appendix C we shall provide further evidence in support of this conjectured relation from 
direct numerical evaluation of the multiple sums. From this relationship we conclude that 
there is only a single independent nontrivial weight 4 function in addition to products of 
Eisenstein series and odd zeta values. 

To order which corresponds to weight 5, the contributions which do not com¬ 

pletely reduce to linear combinations of products of Eisenstein series and ^-values correspond 
to 1 , 1,1 = C's,!,!) -^ 3 , 1,1 = <^ 2 , 1 , 1 ,!, -D 5 = Cl, 1 , 1 , 1 , 1 , and £* 2 , 2 , 1 - Thus, all contributions at 
weight 5 except for £* 2 , 2,1 are expressible as C-functions. By the same type of arguments 
as were used in the case of weight 4, we will be led to conjecture the following algebraic 
relations at weight 5, 

4 OC 2 , 1 , 1,1 = 300C3,i,i + 120E2£3 - 276E5 + 7C(5) 

Ci,i,i,i,i = 60C3 ,i,i + 10E2£3-48E5 + +10C(3)£2 + 16C(5) 

10^2,2,1 = 20C3,i,i-4E5 + 3C(5) (1.11) 

Remarkably, these three relationships show that the four non-trivial weight 5 modular func¬ 
tions may be expressed as a sum of a single non-trivial function in the Ca,h,c series and 
products of Eisenstein series Eg and C-values of odd argument. 

This raises the interesting mathematical question of how many independent C- and D- 
functions there are at any given weight, and what the algebraic and differential structure is 
of the ring of such modular functions. This question is the modular generalization of the 
problem of identifying a basis of independent generators of the polynomial ring of multi-zeta 
values that arise as coefficients in tree-level amplitudes. However, in order to resolve this 
issue it will be necessary to study the expansion of the iV-graviton genus-one amplitude for 
all iV > 4 (the fact that extra diagrams appear at higher weight is clear from the N = 5 and 
N = 6 cases discussed in [39]). This will be discussed further in the hnal section. 

As a consequence of this analysis, we will show that the integrands on moduli space for 
the low energy expansion of the four-graviton amplitude in the Type II string to order [a')^ 
with w < 5, may all be recast in the following general form, 

A{f^) + d^{Eg>,as')) ( 1 . 12 ) 

where A is the Laplace-Beltrami operator in r, is a modular function of weight w, and 
is a polynomial of degree 2 and weight w in the integer weight Eisenstein series Eg' and zeta- 
values C{s'). This form of the integrand allows for an analytical evaluation of the integration 
over r, and thus of the contributions to the effective higher derivative interactions. 

1.2 Organization 

The remainder of this paper is organized as follows. The structure and general features of 
the low energy expansion of the genus-one amplitude are reviewed in section 2. The coeffi¬ 
cients of the higher derivative terms are described in terms of world-sheet Feynman diagrams 
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which reduce to multiple sums over the discrete word-sheet momenta for an arbitrary fixed 
value of r. In section 3, we present a detailed discussion and derivation of the differential 
and algebraic relations obeyed by the modular functions Ca,b,c of arbitrary weight, and we 
exhibit these relations explicitly up to weight 11. In section 4, we introduce the generating 
function W for the modular functions Ca,b,c, derive a differential relation for the Laplace 
operator on W, express this equation in terms of the SO{2, l)-invariant Laplacian, and use 
these ingredients to prove the structure of the inhomogeneous Laplace-eigenvalue equations 
conjectured in section 3. In sections 5 and 6 we discuss the conjectured relations for contri¬ 
butions at weights 4 and 5 respectively, and we present evidence for these conjectures based 
on the asymptotic expansions near the cusp. In section 7 we review the one-loop results up 
to weight w = 3 (order D^IZ^) that were obtained in [9,37], and then make use of the results 
for weights 4 and 5 to derive the full higher derivative interactions up to order by 

performing the r-integral of the corresponding modular functions. In section 8 we end with 
a discussion and comments on higher point functions, higher genus expansions, and toroidal 
compactification. In appendices A and B we review and fnrther develop the calculations 
of asymptotic expansions of non-holomorphic modnlar fnnctions near the cusp. Finally, in 
appendix C we have inclnded a description of a preliminary attempt at a nnmerical check of 
our conjectured relations between the modnlar fnnctions. 


2 Structure of the low energy expansion 

We will be particularly interested in the low energy expansion of terms in the pertnrbative 
expansion of the four-graviton amplitude in either of the ten-dimensional Type II theories in 
flat space-time. In this section, we shall review and clarify both expansions, and give their 
formnlation in terms of world-sheet Feynman diagrams. 

2.1 Perturbative expansion of the four-graviton amplitude 

In either Type II theory, the pertnrbative expansion is in terms of the string conpling gg. In 
Type IIA is the only (real) modulus, while in Type IIB gg is the inverse of the imaginary 
part of the single complex modulus that is sometimes known as the “axion-dilaton”. Type 
IIB is invariant nnder the action of the non-pertnrbative dnality group SL{2, Z) on this 
complex modulus. In either Type II theory, the pertnrbative expansion is given as follows, 

OD 

ki, gg) = ^ A;*) + C>(e"^/®") (2.1) 

h=0 

The exponentially snppressed terms are generated by non-perturbative effects. The h-loop 
amplitnde ki) results from the conformal field theory of matter and ghosts on genus- 

h world-sheet (snper) Riemann snrfaces, and depends on the space-time momenta ki and 
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polarization tensors e* which characterize the four external graviton states labelled by i = 
1,... ,4. Under the Lorentz group S'0(l,9), the ki transform in the vector representation, 
while the e* is a rank two tensor. Total momentum is conserved, the momenta ki are massless, 
and the polarization tensors are transverse. 

As a consequence of space-time supersymmetry of Type II superstrings, the dependence of 
the four-graviton amplitude ki, Qs) on the polarization tensors e* is through a Lorentz 

scalar multiplicative factor, schematically denoted V}. Here, TZ^pucr ~ k^k^ep^ stands for 
the linearized Riemann tensor, and is the unique scalar contraction of four powers of the 
Riemann tensor compatible with maximal supersymmetry.® The remaining dependence of 
the amplitude is through the dimensionless and Lorentz invariant variables Sij = —a'ki ■kj/2 
for 1 < i < j < 4. In view of the zero mass condition for the graviton and total mo¬ 
mentum conservation, the variables Sij obey a number of relations which are conventionally 
parametrized as follows, s = Si2 = <§34, t = S13 = S24, and u = Sn = S23 subject to the 
condition s + t + u = 0 . 

In terms of 7^^, the tree-level amplitude is given by. 


_ kIo n* r(i - s)r(i - t)r(i - u) 
stu r(i-i-s)r(i + t)r(i-i- m) 


( 2 . 2 ) 


where k^q is Newton’s constant in 10 dimensions. 

The genus-one amplitude reduces to an integral over the moduli space Adi of genus-one 
Riemann surfaces of a Lorentz scalar function Bi of the Lorentz-invariant variables s, t, u, 


A^i\ei, ki) = 2 tt / d^i Bi{s,t,u\T) (2.3) 

Jmi 

As usual, we identify A4i with a fundamental domain for the modular group S'L(2,Z) in 
the complex upper half plane, parametrized by the modulus t = ti + iT 2 with ri,r 2 G M, 
0<T2, |ri|<|, l<|r|, and dfxi = dri A (ir 2 /r|. With this normalization the volume of the 
fundamental domain is d/xi = |. Although the variables s,t,u are not independent, we 
shall keep all three as arguments in Bi since Bose symmetry requires Bi to be a symmetric 
function of s, t, u. The partial amplitude Bi is given by the following integral representation,® 


Bi{s, t, u\t) 







(2.4) 


^The explicit form of TZ^ and the expressions for the tree-level and one-loop amplitudes to be given below 
may be found, for example, in [40]. 

®This expression may be generalized to the case of a toroidal compactification of 10-dimensional space- 
time of the form x where T‘^ is a d-dimensional torus specified by a point in the symmetric space 

SO{d, d)/{SOid) (g) SOid)). 
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Here, S denotes the torus with modulus r; the graviton vertex insertion points Zi are inte¬ 
grated over S with volume form d?z = ^dz A dz so that the area of S equals d^z = T 2 ; 
and G{z\t) is the scalar Green function on S given by, 


G{z\t) 


d;(0|r) 



(2.5) 


in terms of the Jacobi 6*-function 9i{z\t). The integral representation for Bi in (2.4) at hxed 
r is absolutely convergent only for limited ranges of s, t, u, but may be analytically continued 
throughout C^, as was shown in [41]. The function Bi can thereby be shown to possess poles 
at negative integer values of s, t and u that arise from singularities in the product of the 
vertex operators for the gravitons. The amplitude that results from the integral over r 
in (2.3) not only has poles in s, t and u corresponding to massive closed string intermediate 
states, but it also has the branch cuts expected by unitarity that arise from the propagation 
of two-particle (massless and massive) intermediate states. 

There are predictions based on supersymmetry and duality for the exact moduli-dependent 
coefficients of the three lowest order terms in the low energy expansion of the four-graviton 
amplitude, namely the 7?.^, and interactions. The explicit expansion of the 

genus-one amplitude up to this order [4, 25] conhrms the predicted values of the one-loop 
contributions to each of these coefficients. 

The genus-two amplitude A^\ei, ki) is also known in terms of an integral representation 
over the moduli space of genus-two Riemann surfaces, along with an integration over the 
four graviton vertex operator insertion points. The coefficient of the hrst term in its low 
energy expansion was evaluated in [8] and shown to match the two loop prediction based 
on space-time supersymmetry and duality for the coefficients of the D^TZ‘^ interaction (an 
analysis using the pure spinor formalism was subsequently given in [42]). Likewise, the 
coefficient of the D^TZ'^ interaction was evaluated in [12] and agrees with its predicted value. 
There is a candidate expression for the low energy limit of the genus-three four-graviton 
amplitude [11], which produces a value for the three-loop coefficient of the TZ'^ interaction 
with an apparent discrepancy of a factor of 3 with the predicted value. 


2.2 Low energy expansion of the four-graviton amplitude 

In brief, the low energy expansion of a string amplitude, such as the four-graviton amplitude 
of A^^\ei,ki,9s), is an expansion in the regime where |sjj| 1 for all i,j, namely when all 
momenta are small compared to the string scale, set by a'. 
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2.2.1 Tree-level 


For tree-level, the expansion may be obtained directly from the explicit formula for the 
amplitude of ( 2 . 2 ), and the expansion of the F-function, and is given by. 




stu 


exp 


E 

, n=l 


2C(2n + l) 
2n -I- 1 


(_g2n+l £2n+l ^2n+l^ 


( 2 . 6 ) 


The low energy expansion of the four-graviton amplitude is seen to be quite special since all 
expansion coefficients will be polynomials in Riemann (^-values with rational coefficients.^ 
Note that only (^-values with odd argument occur in view of the fact that the logarithm of 
the ratio of F-functions in (2.3) is odd under reversing the signs of s,t,u. These amplitudes 
produce higher derivative interactions that are symbolically denoted by where 

stands for the scalar contraction of 2 w covariant derivatives produced by the corresponding 
symmetric and homogeneous polynomial in s, t, u of degree w. A slight rearrangement of the 
expansion is obtained by expressing the symmetric polynomials in s, f, m, in terms of its two 
generators a 2 and 0 - 3 , where. 


a, = s" + r + (2.7) 

subject to the constraint ai = s + t + u = 0. The contribution of weight w then arises as a 
linear combination of monomials crfal whose exponents p, q are related to the weight by, 

w = 2p + 3q (2.8) 

The hrst few terms in this rearranged expansion are as follows, 

Ag \ei, ki) = k\q TZ^ -|- 2d(3) -|- ("(5) 1 T 2 H cg -| ^ CJ 2 -f -C(5)C(3)c’'2cr3 (2.9) 

Because U 2 n+i is divisible by (T 3 = stu for all n > 1, only the leading term in the low energy 
expansion exhibits a pole for small s, f, u. This pole is associated with the exchange of mass¬ 
less states, and may be consistently isolated. The remaining terms contribute respectively 
to the higher derivative interactions of the form, 77"^, D^TZ'^, D^TZ"^, D^TZ"^, and D^^TZ*. 


^The above all-orders low energy expansion of the four-graviton amplitude has been generalized to am¬ 
plitudes with N scattering strings for arbitrary values of N and the resulting coefficients for > 4 are 
generally rational multiples of multi-zeta values of odd weight [28,30]. For weights lower than or equal to 10, 
all multi-zeta values may be reduced to polynomials in ordinary C-values, but above weight 11 irreducible 
multi-zeta values also appear. In the open-string case multi-zeta values enter for weights > 8. This raises 
the interesting question of whether the coefficients in the expansion of the one-loop N = 4-string amplitude 
are again rational multiples of monomials of Riemann zeta values or whether they involve multi-zeta values, 
or more exotic objects. 
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2.2.2 One-loop 

For the genus-one amplitude of (2.3), no poles occur at s, t, or m = 0 thanks 

to space-time supersymmetry, but the amplitude does have branch cuts starting at 0 due 
to the fact that massless strings are being exchanged in closed loops. These non-analytic 
contributions arise from the integration over the genus-one moduli space in (2.3). They 
may be isolated and subtracted to leave analytical terms which provide the higher derivative 
interactions. Thus, a safe point of departure to construct the low energy expansion at one- 
loop is the partial amplitude Bi at hxed modulus r of (2.4). 

The structure of the singularities in s, f, u of the genus-one partial amplitude Bi at fixed r 
is governed by the operator product expansion of the four graviton vertex operators. A 
singularity is produced when two insertion points coincide, namely when Zi — Zj ^ Q for 
i ^ j- The singularities result from local behavior on the surface only, and thus are the same 
for genus-one as they were for tree-level, and consist of simple poles in s, f, u at negative 
integers. Therefore, the integral representation for Bi in (2.4) admits an expansion in powers 
of s, t, u which is absolutely convergent in a region of unit radius in s, t, u. The low energy 
expansion is then simply given by the Taylor expansion of the exponential in (2.4), 


i3i(s, f, u\t) 


£(s,f,M; Zi\T) 


°° 1 1 /• ^ 

/ \\_d^ZiC{s,t,u-Zi\Ty 


SijG{Zi-Zj\T) 


( 2 . 10 ) 


The term of order w contributes to the effective interaction schematically represented by 
Pqj, w < 5, the weight w labels the interaction uniquely, but for weights 

w > 6, degeneracies appear as different Lorentz contractions of the momenta exist, such as 
and (t| for w = 6. These degeneracies are systematically resolved by using a more rehned 
notation for the interaction instead of with w = 2p + 3q. The corresponding 

expansion in powers of <72 and <73 is given as follows, 

CxD 

Bi{s,t,u\T) = al, (2.11) 

p,q=0 

The coefficient functions are modular invariant and have weight w = 2p + 3q, 

following the same notation used for tree-level in (2.8). The basis functions j^^’'^^(t) in this 
expansion follow the conventions in [9,25] and are suited to the four-graviton amplitude®. 


®The more general basis used in [39] generalizes to A^-graviton amplitudes with iV > 4 and facilitates 
comparison of tree-level and loop amplitudes, and hence is adapted to comparison of the genus-one coefficients 
with the multi-zeta values that arise in the tree-level expansion. 
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The coefficient of the interaction erf in the amplitude (2.3) is obtained by integrating 
the modular function over M.i. The resulting expansion of the amplitude therefore 

has the form, 

fc) = + H"-"' tr. + a, + HP-"' tr? + =“’" ffa 

(T2 + fTg + ■ ■ ■ + non — analytic terms ) (2.12) 

The analytic part of this expression resembles that of the tree-level amplitude, with the 
monomials in zeta values in (2.3) replaced by the coefficients 

As indicated earlier, such integrations will in general require regularization and subtrac¬ 
tion of the non-analytic contributions in s, f, u which arise from loops of massless strings. The 
structure and evaluation of such integrals will be discussed in section 7. For the remainder 
of this section we will be concerned with the structure of the modular coefficients j^^’'^^(t), 
and their Feynman diagram representation. Our aim is to determine various properties of 
the functions j*'^’'^^(t). 

2.3 World-sheet Feynman diagrams 

In this subsection we will review the tools needed to study the low energy expansion of the 
genus-one amplitude dehned by (2.3) and (2.4). Each j‘^^’^)(r) function is given by integration 
of the sum of products of tc = 2p -|- 3g Green functions over the toroidal world-sheet with 
complex structure r. This dehnes a sum of Feynman diagrams of weight w, each of which 
is a modular function of r. We will now describe the systematics of this calculation, which 
covers and extends material in [4,25,39], and is a prerequisite for the evaluation of higher 
order interactions in subsequent sections. 

The modular functions j^^’'^)(t) in the expansion of (2.11) are given by (2.10) in terms 
of integrals over four copies of the torus S of products of w scalar Green functions. To 
evaluate these integrals, it will be convenient to use an alternative representation of the 
Green function G which renders its double periodicity fully transparent. We parametrize a 
point 2 : on the torus S of modulus r by real coordinates a, /3 ^ M/Z, 

z = a + Tf3 (2-13) 

The Green function G is then given by the Fourier series over integers m,n E Z, 

G{z\t) = Go{t) + Q{z\t) (2.14) 


with 


g{z\T)= 

(m,n)^(0,0) 


( 2 . 15 ) 
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which may also be viewed as the definition of a Kronecker-Eisenstein series. The Fourier 
modes ^(m, n|r) are given by, 


Q{m, n|r) 


r2 

nlrriT + nP 


(2.16) 


for {m,n) ^ (0,0) and we shall set ^(0,0|r) = 0. The integers m and n may be viewed 
as the components of the world-sheet momentum of the scalar field, which will be denoted 
p = mr -|- n, and we shall sometimes use the alternative abbreviated notation ^(p) = 
Q{m,n\T). The zero mode contribution Go{t) = 2 In (27r |? 7 (r)|^) will drop out of (2.10) 
and any perturbative string calculation in view of the momentum conservation relation 
s + t + u = 0. This has the immediate consequence that the term linear in G will vanish. 
Therefore, the world-sheet momentum space Green function Q contains only non-zero modes. 

Expanding C{s,t,u-, Zilr)^ of (2.10) in powers of the Green functions gives a monomial 
of weight w in which different powers of G{zi — Zjlr) link the four vertices at Zi, Z 2 , z^, 
Z 4 . Substituting the Fourier series (2.14) for each Green function and performing the Zi 
integrations leads to a sum of world-sheet Feynman diagrams in which each propagator 
joining vertices i and j denotes a factor of G{pij) with a sum over all the (discrete) world- 
sheet momenta p^ subject to momentum conservation at each vertex. 


2.4 Nomenclature for world-sheet diagrams 

In the following it will prove convenient to symbolize these Feynman diagrams in different 
fashions, according to the context in which they are discussed. The problem is that there is 
no simple way of labeling the function corresponding to a general (and complicated) diagram 
- the simplest representation is in terms of the diagram itself. 

A systematic way of labeling an arbitrary diagram, at any order in the low energy ex¬ 
pansion, and with any number N of vertices and a total of w propagators is as follows. Each 
vertex is labelled by an integer i,j = 1,2,..., N, and the number of propagators that joins 
the pair of vertices® {i,j) by iij. Each propagator joining a given vertex pair {i,j) is assigned 
a label {i,j, \ij) where 1 < Xij < iij, and carries world-sheet momentum 
In that case the modular function associated with the diagram is given by the product of 
propagators summed over word-sheet momenta, 

H 11 ^ (Ph') 11 fe £ P^“) (2.17) 

I l<i<N \ Xik = l / 

iPij s 

where the product of (5’s enforces conservation of momentum at each vertex. 

Hi iij were restricted to take only the values 0 or 1 it would be the incidence matrix for the graph - 
however, in general we have iij G N. 
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This notation has redundancies since the modular function associated with the diagram 
is invariant under an arbitrary relabeling of the vertices. It is also needlessly complicated 
for many purposes, so, in the following we will label the functions associated with diagrams 
in a manner that is adapted for the particular problem at hand, bearing in mind the general 
expression in (2.17). 

2.5 Diagrammatic representation of modular coefficient functions 

In this section, we will be interested in 

• those diagrams that contribute to the four-graviton amplitude and, 

• diagrams belonging to a particular infinite subset of all the diagrams that produce 
higher-point functions and which we can determine for any weight. 

2.5.1 Contributions to the four-graviton amplitude 

In this case a general diagram has iij propagators joining any pair of vertices labelled i and j, 
where i,j = 1, 2, 3,4. Therefore there are six {i,j) pairs so a general diagram has ii,..., Eq 
propagators. The notation we will adopt is summarized by the following diagram: 



where the label i on the link ^ indicates the number of propagators joining the 
corresponding pair of vertices. The shape of this diagram is a tetrahedron that is symmetric 
in all edges and vertices and its weight is given by the number of propagators, so that 

6 

w = J2^r (2.18) 

r=l 

where 0 < C- When C = 0 for particular values of r, the diagram degenerates to a simpler 
diagram. In such cases, we shall omit the index that vanishes. 

Certain classes of diagrams automatically give vanishing contributions. Thus, any dia¬ 
gram that has a vertex Zi connected to a single propagator is zero due to the fact that the 
integral of Q{z\t) over vanishes (so that there is no zero momentum mode in ^(m, n|r)). 


i°For example '^^ile £>^,o,o,o.;o.o = Dg, etc. 
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Likewise, one-particle reducible diagrams vanish since any such diagram contains a zero- 
momentum propagator. Furthermore, diagrams involving the self-contraction of a propaga¬ 
tor do not arise since such diagrams have a coefficient proportional to = 0, where kr is 
the momentum of an external graviton. 

The following is a selection of the non-vanishing diagrams that contribute to the expansion 
of the four-graviton amplitude. 

(1) Diagrams involving a pair of vertices. 

When all but one of the £r = 0 the results is diagram with £ propagators joining a pair 
of points, represented by 



= n, 


(2) Diagrams involving three vertices. 



(3) A few of the diagrams involving four vertices. 



In calculating the contribution of the diagrams of a given weight to a term of given 
transcendentality in the low energy expansion we need to sum the diagrams of that weight 
with combinatorial coefficients that were evaluated up to weight 6 in appendix C of [25], and 
can easily be evaluated to an arbitrarily high weight. 
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2.5.2 A subclass of diagrams contributing to amplitudes with higher N 

The diagrams we have discussed so far are ones that contribute to the four-graviton ampli¬ 
tude, They are also a sub-class of the diagrams that contribute to iV-graviton amplitudes 
with N > 4. However, when N > 4 there are more diagrams for two distinct reasons: (i) 
There are more vertices so there are N{N — l)/2 ways of joining a pair of vertices; (ii) There 
are 2(iV —4) momentum numerator factors in any diagram. Many of these can be eliminated 
by integration by parts, thereby reducing the diagram to one with fewer (or no) momentum 
factors. This systematics was considered in detail in [39], for the N = 5 and N = 6 cases 
where the diagrams that contribute up to weight 6 were enumerated. 

The evaluation of the complete set of diagrams for the N > 4 amplitudes will not be 
pursued here. However, in section 3 we will consider an inhnite subset of diagrams, in which 
single propagators are joined to form chains of the form 





a 


A single such chain is zero because of the rule that a diagram with a vertex connected to 
a single propagator vanishes. The hrst nontrivial case is that of two chains of length a and 
b, which are represented by. 



This is equivalent to a circular closed chain of length a + b. In an obvious generalization 
of the earlier notation this could be denoted i('r), but this is clumsy and does not 

a + h 

generalize to cases with more chains. For that reason we will use a different notation, in 
which a graph with a number of chains with lengths ai,... ,ap is denoted Ca^^...^ap(r). From 
the above hgure it easily follows that. 


CUr) = E 

(m,n)^(0,0) 


T-2 


TT Imr -|- nP 


a+6 




(2.19) 


which is the Eisenstein series of weight a + b. 

Signihcant and very interesting challenges arise in determining the functions associated 
with three or more chains. The case we will concentrate on in section 3 will be graphs with 
three chains, which are denoted Ca^b,c with a graphical representation. 
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In the next section we will derive a set of coupled inhomogeneous Laplace-eigenvalue equa¬ 
tions that these coefficients satisfy, which will lead to solutions for arbitrary a, b and c. 

There is a variety of obvious generalizations of these classes of diagrams, such as the 
functions Cai,a 2 ,a 3 ,...,ap already dehned in the introduction in (1.1), that has p chains of 
arbitrary length. Another obvious generalization is to replace the single propagator in each 
link of a chain by a self-energy bubble, so that a single chain has the form. 



This self-energy bubble diagram will be a key ingredient in the calculation of the diagram 
^ 2 , 2,1 by direct summation in appendix B. 


2.6 The Laplacian on modular coefficients via deformation theory 


As explained in the Introduction, it is the Laplacian A on the genus-one moduli space Ati 
acting on the partial amplitude Bi in (2.10) that will be of central interest in this paper. 
In terms of the standard parametrization of At i by a modulus r = Ti + ir 2 in the complex 
upper half plane, the Laplacian takes the form A = and therefore involves a mixed 

holomorphic - anti-holomorphic derivative on moduli. Such derivatives can be evaluated 
with the help of the theory of deformations of complex structures in terms of associated 
Beltrami differentials fi and Ji which, for the torus E, may be chosen to be constant on E. 

The Beltrami differential parametrizes an inhnitesimal deformation of the flat metric \dz\‘^ 
on a torus with modulus r to a metric \dz + p-dzl"^. Under the quasi-conformal transformation 
z' = {z + p)/{I + p) the torus with modulus r and metric \dz + pdz\^ is mapped to a torus 
with modulus r' = {t + /if)/(l -|- p) and metric \dz'\^ (up to a constant Weyl factor which is 
immaterial thanks to conformal invariance of the amplitudes). The Laplacian A may then 
be simply expressed in terms of deformations by p via the formula. 


A = d^d-^ (2.20) 

In this formalism, it is the metric that is being deformed, but the coordinates that are left 
unchanged. Thus, the only thing that needs to be varied in the modular function Bi of (2.10) 
is the Green function Q. Its deformations to hrst order in p and p are given as follows, 


d^,Q{zi 



— / d^zd^Q{z - Zi)d^G{z - Zj) 

TT ./v 


( 2 . 21 ) 


this section, to shorten the notations, we shall exhibit neither the dependence of Q^Bi and C on r, 
nor the dependence of Bi and C on s, t, u. 
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along with the complex conjugate equation. The second order mixed variation d^dfiQ{zi — Zj) 
evaluated at /i = 0 vanishes in the interior of moduli space. 

Therefore, it is now straightforward to express the Laplacian of Bi in terms of these 
deformations, and we hnd. 




w\ 


To 

w=2 ^ ^ i=l 

The deformation formula for £ is given by. 


]Jdh,£--^(d,£)(dp£) 


( 2 . 22 ) 


1 ^ /I ^ 2 -j 


(fzd^Q{z - Zi)dzQ{z - Zj) 


(2.23) 




At low weight orders, this formula may be simplihed and worked out quite explicitly. We 
consider, for example, the case of the Feynman diagrams given by 


Dp = 


(Pzi f (Pz2 


Q{Zi - Z 2 Y 


(2.24) 


Jt. T2 Jt. a 

Applying these techniques yields their Laplacian, and after some simplihcations, we obtain, 

4 


AD, = + 


n / “ Z-i)Q{z2 - Zi) 

^ D is 

xdzsG{z2 - zY)dz^Q{zi - Z4)d^^d2^Q{zi - Z 2 ) 


(2.25) 


?-2 


When i = 2, the integral on the right hand side vanishes in view of the last factor of the 
integrand, and we obtain AD 2 = 2 D 2 in accord with the fact that we already know D 2 to be 
simply equal to the Eisenstein series E 2 . For i = 3, the first term vanishes, and the integral 
may be simplihed since the last factor of the integrand evaluates to 


TT 


< 9^1 < 952 ^( 2 : 1 , Z2) = ttS{zi - Z2) - 

T 2 


(2.26) 


The integration over the (5-function term in (2.26) sets zi = Z 2 so that the integrand becomes 
a total derivative in Z 3 and Z 4 and the integral of this contribution vanishes. The integration 
over the constant term in (2.26) gives. 


6 ^ f 

AD 3 = —3 TT / d^Zi G(zi - Z 3 )g(z 2 - Z 4 )d^^Q{z 2 - Z 3 )d 2 ^Q{zi - Z 4 ) 

ZA is 


(2.27) 


Transforming the derivative in .24 into a derivative in zi by translation invariance, and then 
integrating by parts in both zi and ;23 exposes the combination dz^^dz^Gizi — Z 3 ) which is 
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evaluated with the help of (2.26). This time the contribution from the constant term cancels 
out, and its is the ^-function which survives and gives, 

AT >3 = — TT / g{z 2 - Zs)g{zs - Z 4 )g{z 4 - Z 2 ) (2.28) 

^2 i=2d^ 

The right hand side readily evaluates to an Eisenstein series, and we find ADa = 6 E 3 , which 
will also be arrived at by a slightly different procedure in the next section. 

The evaluation of higher order contributions by this method has not yet been system¬ 
atically investigated, as we shall use more algebraic techniques to go to higher order in the 
subsequent sections. 


3 The modular functions C, 


a,b,c 


This section will be devoted entirely to the study of the modular functions Ca,b,c which 
corresponds to the case p = 3 in (1.1). The functions Ca^b,c evaluate an inhnite class of 
Feynman diagrams with three chains, as explained in section 2.5. It is for this case that 
we will be able to obtain the most general and the most precise results, for which we shall 
present complete proofs in this section, and the subsequent section 4 with the help of a 
generating function for Ca,b,c- 

To start, we consider a generalization of the set of functions Ca,b,c dehned as follows, 

^ dnifldn,oimin2 - nim2Y^ (r2/7r) 


CibA 


T = 


“ |mir nip“|m2r n2p^|m3r 

jTl'f jy—U 


(3.1) 


for any non-negative integer k. The original function Ca,b,c corresponds to the special case 
k = 0. The Kronecker 6 factors enforce the vanishing of the total momentum components 
m = mi + m 2 + m 3 and n = ni -|- n 2 + ns. Of interest here will be the case where a, b, c are 
integers such that the sums are absolutely convergent, but the functions dehned 

for complex values of a, b, c as well. The functions Ca^b,c are invariant under permutations of 
the indices a, b, c for all k in view of the fact that momentum conservation m = n = 0 implies 
{min 2 — nim 2 )^ = {m 2 ns — n 2 msY = {jnsUi —nsmi)"^. In the following we will only consider 
integer values of a, b and c, and for the most part we will follow the convention of using this 
permutation symmetry to label the indices so that a>b>c. Finally, the functions 
are modular invariant because the combination {min 2 — nim 2 Y is a symplectic invariant. 


3.1 Algebraic representation of the Laplacian 

The key tool that we shall use in the analysis of the modular functions i® f^ct that 
the Laplace-Beltrami operator A on the complex upper half plane, 

/^ = rUdl+dl) (3.2) 
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acting on the functions admits an algebraic representation. Indeed, for a function 
of weight tc = a + 6 + c, we hnd the following result by explicit calculation, 

(a - - 4fcc‘+j,,,+, - ( 3 . 3 ) 

The functions for successive values of k are related by the following algebraic identity, 

A^k-\-l _ I ^k I ^k A\ 

~ ^a-2,b,c ^a,b-2,c ^a,b,c-2 ^^a,&-l,c-l ^^a-l,b,c-l ^^a-l,b-l,c 

which is easily proven by applying the relations for momentum conservation m = n = 0 to 
derive the following polynomial identity valid under the summation in (3.1), 

3 

\mrT + rirl"^ — \mrT + nr\‘^\mr'T + = —4r2 (min2 — nim2f (3.5) 

r=l r^r' 

The results of (3.3) and (3.4) may be combined to eliminate the functions with upper index 
/c + 1 , and produce a relations between functions of upper index k only, 

(a - o(o - 1) - 6(6 - 1) - c(c - l))c*t., 

— +<i6^cy^ + cy,^ + cy, — 2cyj ,^- 1 ) 

+'>c(cy2,^+i,2+i+- 2cyi,j,,+, - 2cy, 

+ca(c*_,,,,2+i + cy,,,_2,,+, + - 2 C*j_i,2+i - (3,6) 

This expression is closed on the functions C^.. and we may set A; = 0 to recover the functions 
in which we are interested. The application of the Laplacian is now reduced to an algebraic 
operation. Note that, although the functions of central interest in string theory have positive 
integer a, b, c, the above formulas are valid for a, 6 , c G C as long as their values are such that 
all functions involved correspond to convergent sums. 


3.2 Evaluating Ca,6,o and Ca,b-i 

The right hand side of (3.6) may end up involving a lower index which vanishes or which 
equals -1. Both of these functions will appear in applying the Laplace operator to functions 
Ca,b,c with a,b,c > 1. We shall now evaluate these terms separately. The corresponding 
sums then simplify considerably, and reduce to. 


C«.ft.n(r) = 


CaA-iir) = 


(mr,nr)^(0,0) 


(mr,nr)y^(0,0) 


7r“+”|mir + nip“|m2T + n2f“ 

I (mi + m 2 )r + ni + n 2 | 

7r“+^“^|mir + niP“|m2T + n2h 
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In Ca,b,o, the sum excludes the values (m 2 ,n 2 ) = (—mi,—ni); restoring their contribution 
allows evaluating the sums in terms of Eisenstein series. In Ca^b,-i the exclusion of the values 
(m 2 ,n 2 ) = (—mi, —ni) is automatic thanks to the presence of the numerator, and decouples 
the sums into a product of Eisenstein series. The final results are as follows, 

C'a,b,0 — Ea Eb Ea-\-h 

Ca,b-1 = Ea-l Eb + Ea Eb-l (3-8) 

Having established these special values we will now proceed to determining the equations 
satisfied by the modular functions Ca,h,c- To start with, we will study the functions with 
weights w = 3,4, 5, before analyzing the case of general weight w. 


3.3 Laplace equation for Ca,b,c at low weights 

Of central interest in string theory applications will be the functions Ca,b,c with integer 
a,b,c> 1. We shall now derive explicitly the Laplace operator on the corresponding functions 
for tc = a + 6 + c = 3,4, and 5. 


3.3.1 The case Oi 11 = D 3 

In this case we have a = & = c = 1. Formally, the Laplace equation derived from (3.6) 
becomes AOi^i^i = 302,2,-i — 602,i,o- Actually, both functions on the right hand side of this 
equation are divergent, which may be attributed to a logarithmic divergence of the sum over 
world-sheet momenta, as may be seen by evaluating them with the help of (3.8) and the fact 
that El is divergent. 

To carry out the evaluation, we combine the sums of both terms,and observe that the 
combined sum is convergent, and given by. 


^*2,2,-! ~ 2(72,1,0 — 


{mr,nr)^0 


((miT -F ni)(m2r -F n-2) + c.c.) 
TT^ ImiT -|- ni|"^|m2r -|- n2\^ 


(3.9) 


where we still use the notation m = mi -|- m2 -l- m3 and n = ni + n2 + n^- The sum 
excludes the terms with (mi,ni) = (—m2, — ^2). Restoring this contribution to the double 
sum makes each factor in the double sum odd under {mr,nr) —)■ {—nir, —fir). Hence the 
unrestricted double sum vanishes, leaving only the contribution from (m 3 , n^) = ( 0 , 0 ), which 
gives ( 72 , 2 ,-! — 2 ( 72 , 1,0 = 2 E 3 . As a result, the Laplace equation becomes. 


A(7i,i,i — 6E3 


(3.10) 


alternative procedure is to regularize all functions in the equation, for example with the help of a 
cutoff A on the world-sheet momenta |m|,|n| < A. The divergent function Ei then gets replaced by its 
momentum cutoff version which will cancel out of all calculations of the Laplacian. An alternative 
regularization is by continuing the functions Ca,b,c to non-integer values by infinitesimal shifts a ^ a + e etc. 
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Using ( 1 . 4 ), we recast the equation as A — E^) = 0 . Its solution is given by, 

Ci,i,i = .£^3 + C( 3 ) ( 3 . 11 ) 

where we have used knowledge of the power-behaved terms in the large-r2 expansion of Cm, 
given in equation (B.8) in [ 25 ], 

Cu,i = ^ + C(3) + Ah + (3,12) 

In this equation, and many of the following, we have defined the quantity, 

y = 7lT2 ( 3 . 13 ) 

which eliminates many factors of tt from subsequent formulas. The result ( 3 . 11 ) has been 
obtained earlier by Don Zagier [ 38 ] by direct calculation of the sums over m, n. 

3.3.2 The case ^ 2 , 1,1 = D 244 

Using the expression for the Laplacian in ( 3 . 6 ), we readily hnd again a formal relation, 

(A — 2 ) 6*2,1,1 = —+ 46 * 3 ,2,-1 — 56 * 2 , 2,0 ( 3 . 14 ) 

where the expressions 6*3,1,0 and 6*3,2,-! are separately logarithmically divergent through 
terms linear in Ei, which require regularization to E^. In the combined contribution we 
may directly apply ( 3 . 8 ), resulting in the cancellation of all dependence on E^. The Laplace 
equation then reduces to, 

(A - 2 ) 6*2,1,! = 9^4 - E^ ( 3 . 15 ) 

This equation has the same structure as the inhomogeneous Laplace-eigenvalue equation 
that arose in a rather different context in [ 9 ]. 

3.3.3 The cases 62 , 2,1 “ ^i,!,!,!;! 6311 = 6 ) 2 , 1 ,!,! 

We will now consider the two weight -5 modular functions 62,2,1 and 63,1,1. Using the expres¬ 
sion for the Laplacian ( 3 . 6 ), we readily hnd the relations, 

A62,2,1 = 463,3,-1 — 863,2,0 

(A — 6)63,1,1 = 362,2,1 — 664,1,0 ~ 1063,2,0 + 664,2,-1 ( 3 . 16 ) 

Again, the relation on the second line is formal and requires regularization as in the earlier 
cases. Substituting the corresponding expressions from ( 3 . 8 ) gives, 

A62,2,1 = 865 

(A — 6) 63,1,1 = 362,2,1 + I66/5 — 46/26/3 ( 3 . 17 ) 
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(3.18) 


In particular, using /S.E^ = 20i?5, it follows that A(C2,2,i — f-^s) = 0, and thus, 


C*2,2,l — 
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where the coefficient of the C(5) term is determined from knowledge of the terms that are 
power behaved in the limit y ^ oo, given in [25]. Substituting this in (3.17) leads to. 


(A - 6) ^3,1,1 = y^5 - 4E2E3 + ^ 
which is analogous to (3.15) and will be discussed further in section 6. 


(3.19) 


3.4 Poincare series representations 

The inhomogeneous Laplace eigenvalue equations satished by the modular functions dis¬ 
cussed in this paper have a structure analogous to the equation that arises in the context of 
S-dual properties of the coefficient of the interaction obtained in [9]. The solution to 

this equation was obained in the form of a Poincare series in appendices A and B of [10]. 
Here we will derive a similar Poincare series representation for the specihc example of the 
modular function 6*2,1,1, but the same methods can be applied to obtain solutions for any of 
the Laplace equations that arise in subsequent sections. 

We begin by recalling the classic Poincare series representation for the Eisenstein series, 
E,(r) = (l'^7(^)) (3.20) 

7eroo\5L(2,z) 

where 7 acts on r by Mobius transformations, 

= ^ 7=(: ^,)€SL(2,Z) (3.21) 

and Poo the subgroup of SL{2, Z) which stabilizes r = ioo. Note that we can also represent 
any of the modular functions Ca,b,c by a sum over 2x2 matrices, 

= E (detT)“-“-''-'(lm7(r))”(lmS7(T))‘(lmST7(r))' (3.22) 

7eM(2,Z) 

where S{t) = —l/r, T(t) = t + 1 and M(2, Z) is the set of 2 x 2 matrices with integer 
entries, and non-vanishing determinant. Although this expression is not a Poincare series it 
can be recast as one in the following manner. 
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We will also use the Fourier series for a non-holomorphic Eisenstein series, valid for any s, 


E.{t) 


2((2s) , 2r(i, - l)C(2s - 1) 

-r* r(i,K-4 ^ 

AV.(27r|Af|T2)e"‘’"" 

AT^O 


(3.23) 


where Ky{z) is the modihed Bessel function of order z/, and the divisor sum is given by 
<^k{n) = 'Ylid\n positive divisors d. With our normalisation, the Eisenstein series satisfies 
the functional equation T{s)Es{t) = r(l — s)Ei_s(r) for all r in the upper half-plane, and 
s G C. For integer values n of s, the function Kj^_i reduces to a spherical Bessel function, 
and we have the following simplified form, expressed in terms of g = and y = 7 rr 2 , 


En{r) = 




( 2 n)! 


(n — 2 )!(n — 1 )! 


(3.24) 


n—1 


+ 


(n — 1 )! 




N=l 


—^ m\ in — m — 1)\ 

m=0 ' 


where i? 2 n are Bernoulli numbers. We see that when s = n is an integer the series in m 
in (3.24) truncates to a hnite sum with n terms and 0 < m < n for any value of N. The 
purely power behaved terms, together with the leading exponential corrections (the terms in 
the sum with iV = 1) will play an important part of our considerations in sections 5 and 6 . 

Consider now the inhomogeneous Laplace equation for 6 * 2 ,ly in (3.15), which we will 
rewrite as 

E, 


Ea - E^ 


(3.25) 


(A — 2) ( 6*2,1,1 


where we have subtracted 2 OE 4/3 from each side (using the fact that (A — 2)^4 = IOE 4 ) 
so that the source term on the right hand side now has leading behavior 0{y) in the limit 
y —)■ 00 , as is seen from the expansion (3.24) for the cases n = 2 and n = 4. 

We shall now hnd a solution to (3.25) expressed as a Poincare series of the form 

C' 2 ,i,i(r) - ^E 4 (r) = $( 7 r) (3.26) 

7eroo\5L(2,z) 


To determine the function ‘h(r) we recall the Poincare series representation of the Eisenstein 
series of (3.20). Substituting this expression for the cases s = 2 and s = 4 in the source term 
on the right-hand side of (3.25) and substituting the expression (3.26) on the left-hand side 
of (3.25) leads to the inhomogeneous Laplace equation for $ 

4 2 

(A-2)$(r) = ^- —E 2 (t) (3.27) 

^ ^ ^ ^ 2025 45 ^ ^ ^ ^ 
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where the leading behavior of the right hand side is of order 0{y) as y —)■ cx). 

Noting that and g™ are annihilated by the Laplace operator, Ag™ = = 0, we 

seek a solution to (3.27) of the form 

CXD 

$(r) = $o(l/) + $n(|/) (g" + r) (3.28) 

n=l 

Given that lim^^oo *h(r) = 0{y), it follows from Lemma 2.9 in [10] that *h(r) must satisfy the 
boundary conditions lim^^^o‘^n(l/) = 0{y^) for all n > 0. This determines that the solution 
of (3.28) is uniquely given by 



and therefore the Poincare series representation of 6 * 2 ,i,i is determined by substituting <l>( 7 r) 
in (3.26). 


3.5 Laplace equations for Ca^b.c at any weight 

For general weight tc = a + 6 + c>5 with a,b,c > 1 the Laplace equations of (3.6) involve 
several functions Ca,b,c- The corresponding system of coupled equations can be diagonalized 
in order to determine the spectrum of eigenvalues and eigenvectors of the Laplacian acting on 
the space of Ca^^c for weight w. With the help of MAPLE we have carried out this analysis 
up to weight w = 34 included. The results may be summarized by the following, 


3.5.1 Theorem 1 : Structure of the Spectrum 

(a) The eigenvalues of the Laplacian acting on the space of functions Ca,b,c with weight 
tc = a + 6 + c for a, 6 , c > 1 are given by s{s — 1 ) where s = w — 2 m and m runs over all 
integers in the range. 


1 < m < 



(3.30) 


(b) The multiplicity Ug of the eigenvalue s(s — 1) equals 

Us = (3.31) 


Here, [x] denotes the integral part of x, dehned to be n G Z for x in the range n < x < n+1- 
A complete analytical proof will be given in section 4 with the help of a generating function 
W for the modular functions Ca,b,c- 



3.5.2 Illustration of the spectrum up to weight w = 12 

To illustrate the above result on the structure of the spectrum, we shall list below the allowed 
values of s for weights w in the range 3 < tc < 12, along with the eigenvalue s(s — 1) and 
the degeneracy Ug in the notation s(s — 


w = 3 

s = 1 

0 (1) 

w = 4 

s = 2 

2 (1) 

w = 5 

s = 1,3 


w = 6 

s = 2,4 

2 (1) 0 12(2) 

w = 7 

s = 1,3,5 

0 ^^^ © 6^^) © 20^2) 

w = 8 

s = 2,4,6 

2 (^) © 12(2) ^ 3 q(2 ) 

w = 9 

s = 1,3,5,7 

0 (^( ©6(^( ©20(2) ©42(^( 

w = 10 

s = 2,4, 6, 8 

2 (^( © 12(2) ^ 3 q(2 ) ^ 3g(3) 

tc = 11 

s = 1,3, 5, 7, 9 

0 (^( © 6(^( © 20(2( © 42(^( © 72(^( 

w = 12 

s = 2,4,6,8,10 

2(^( © 12(2) ^ 3 q(2 ) ^ 3g(3) ^ gQ(4) 


MAPLE also provides the full inhomogeneous Laplace-eigenvalue equations. We shall list 
these in the subsections below up to eigenvalue 12, which is the lowest eigenvalue at which 
the multiplicity exceeds one. 

3.5.3 The full inhomogeneous eigenvalue equations for eigenvalue 0 

For eigenvalue 0, namely s = 1, the multiplicity is one for each odd weight (and zero for each 
even weight), and we list the corresponding inhomogeneous Laplace-eigenvalue equations 
below for each odd weight, up to weight tc = 9, 



= 6E3 


AC'2,2,1 

= SEs 


^ (^*3,3,1 + (^3,2,2) 

= I8E7 


A (9(74^4^1 + 18(7*4^3^2 + 4(73^3^3) 

= 288Eq 

(3.33) 


The first two cases are already familiar from our analysis at low weights in (3.10) and (3.17). 
These equations were integrated with the help of the eigenvalue equation = s{s — l)Es 
for the Eisenstein series and the integration constant was hxed by matching the asymptotics 
at the cusp. 

Clearly, the equation at odd weight w associated with eigenvalue 0 exhibits a right hand 
side proportional to the Eisenstein series and may thus be integrated to obtain a purely 
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algebraic relation between the Ca,b,c at weight w = a + h + c. The algebraic equations 
corresponding to the differential equations of (3.33) are given as follows, 

Cl,1,1 = C3 + (C(3) 

30C*2,2,i = I 2 T /5 + C(5) 

14(73,34 + 14(73,2,2 = 3i?7 + gi 

9(74,4,1 + 18(74,3,2 + 4(73,3,3 = 4i?9 + gfg (3.34) 

where gj, gg are integration constants which are determined by matching the asymptotic 
behavior of both sides of the equation at the cusp. A general formula for these algebraic 
relations valid for any odd weight will be derived in section 3.9. 


3.5.4 The full inhomogeneous eigenvalue equations for eigenvalues 2 and 6 

For eigenvalue 2, which has multiplicity 1 for every even weight starting at weight 4, we have 

(A — 2)(72,i,i = 9 E 4 — E 2 
(A — 2) (4(73,2,1 + ( 72 , 2 , 2 ) = 52Eq — AEI 
(A — 2) (6(74,3,1 + 3(74,2,2 + 5 ( 73 , 3 , 2 ) = 153Eg — 9E‘^ 

(A - 2) (8C'5,4,i + 8C'5,3,2 + 7C'4,4,2 + 8C'4,3,3) = 336^10 - 16E| (3.35) 

The first equation is already familiar from (3.15). For eigenvalue 6, which has multiplicity 1 
for every odd weight starting at weight 5, we have, 

(A — 6 ) (2(73,1,1 + (72,2,i) = 407^5 — 8 E 2 E 3 
(A — 6) (6(74,2,1 + 2(73,3,1 + 3 ( 73 , 2 , 2 ) = IbQE-j — 2 AE 3 E 4 
(A — 6 ) (8(75,3,1 + 4(75,2,2 + 3(74,4,1 + 10(74,3,2 + 2 ( 73 , 3 , 3 ) = 384F^9 — A 8 E 4 E^ (3.36) 


3.5.5 The full inhomogeneous eigenvalue equations for eigenvalue 12 


For eigenvalue 12 , the multiplicity is 2 for every even weight starting at weight 6 . At weight 


6 and 8 we have respectively. 


(A - 12) 

6(73,2,1 — (72,2,2 
6(74,1,1 + ^2,2,2 _ 

= 

108^6 - 36.F| 
120^6 + El - 36E2.F4 

(A - 12) 

2(74,3,1 + (74,2,2 
8(75,2,1 + (74,2,2 + 2(73,3,2 

= 

66^8 - l 8 El 
190^8 + 18El - A 8 E 3 E 

and at weight 10, 





(A-: 

1 2'j 12(75,4,1 + 12(75,3,2 + 3(74,4,2 + 2(74,3,3 

40(76,3,1 + 20(70,2,2 + 16(75,3,2 + 9(74,4,2 + 14(74,3,3 _ 


(3.37) 


624^10 - 144^2 
1656^10 + 144E| - 360^4.^6 


(3.38) 
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3.6 Structure of the eigenspaces of the Laplacian 

In this section, we shall infer the strnctnre of the eigenspaces of fnnctions Ca,b,c from MAPLE 
calculations up to weight w = 34. Recall from Theorem 1 that the eigenvalues of the 
Laplacian are given by s(s — 1), where s = w — 2m and the allowed values for m are given 
by m = 1, 2 ,..., [(tc — l)/2], and corresponding eigenspace V(tc, s) has dimension, 


dim V(tc, s) 


s + 2 
3 


(3.39) 


In addition to this general structure of the spectrum, one may provide a more detailed 
description of the eigenfunctions leading to statements that will be proven in full using the 
generating function in section 4. 


3.6.1 Structure of the space V(tc, 1) for eigenvalue 0 

For zero eigenvalue, we have s = 1 and therefore the weight is given by tc = 2/r + 1 with 
m = /i for any positive integer fx. The corresponding eigenspace of eigenvalue zero is then 
generated by the following modular functions. 


C/i—p,/i—p—q,2p+q+l 0 < p, 0 < q. 

The range may be explicitly parametrized by 


p = 0,1 ,..., 


/r 


q = 0,1,..., 


3p + 2q < yu — 1 


/i — 1 — 3p 


(3.40) 


(3.41) 


Clearly the structure depends on m modulo 6. 


3.6.2 Structure of the space V(tc,s) for arbitrary w,s 

For odd weight tc = 2/r + 1 > 3 we have odd values s = 2i + 1 with i < fj, and xn = fx — i. 
The eigenspace is restricted by the condition that only the functions Ca,b,c can contribute 
for which the value of a does not exceed fx + i. This set may be parametrized as follows, 

p—q,2p+q+l—2£ 0 ^ P) 0 < q 

3p + 2q < /i + 2£ - 1 (3.42) 

For even weight w = 2fx > 4 and thus even values s = 2i with i < fx, the eigenspace is again 
restricted by the condition that only the functions Ca,b,c can contribute for which the value 
of a does not exceed fx + i. This set may be parametrized as follows, 

p—q,2p+q—2£ 0 < p, 0 < q 

3p + 2q < yu + 2£ 
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(3.43) 



3.7 Embedding eigenspaces at successive weights 

We shall establish the existence of a map, denoted which embeds an eigenspace V(tc,<s) 
into the eigenspace V{w + 2, s) of the same eigenvalue but weight increased by 2, 

C: V(w,s) ^ V(w + 2,s) (3.44) 

for all s = tc — 2m with m = 1, 2,, [(m — l)/2]. The map on V(m, s) may actually be 
dehned on the larger space of all functions with a,b,c> 1 and k > 0, via the relation, 

e {Cl,Jr)) = abCl,^,^,jT) + (3,45) 

The key observation in the proof of the above assertion is that the functions satisfy 

the same Laplace-eigenvalue equation (3.6) for all value of /c > 0. To see this, let a linear 
combination Cj with coefficients J^a,b,c of ^(r)-functions at hxed weight w = a + b + c 
satisfy the eigenvalue equation for one value of /c, 

(a - i,(i, - 1 ))cJ(t) k 0 (t) = J.A.CjJr) (3.46) 

[a,b,c] 

Here, ~ 0 stands omitting inhomogeneous terms linear and bilinear in Eisenstein series, and 
the sum in [a, b, c] is over all a > 6 > c > 1 with a + b + c = w. If (3.46) is satished for one 
value of /c > 0, then it will be satished for all values of k with the same coethcients Ja,b,c- 
Now on the one hand, Cj has eigenvalue s(s — 1) as written in (3.46), while on the other 
hand from (3.3), we have, 

(A-«,(«,-l))CyT) = -4£(Cy‘) (3.47) 

Applying {A — w{w — 1)) to the left hand side of (3.46), and using the fact that this operator 
commutes with the operator (A —s(s —1)), we obtain the relation (A —s(s —1))(£(C^’^^) ~ 0. 
Since the same Laplace-eigenvalue equation is obeyed irrespective of k, it follows that we 
also have the eigenvalue equation, 

(A-s(s-l))e(C'^) ^0 (3.48) 

where it is understood that k = 0, and the weight of ^{Cj) is m -|- 2. 

3.8 Corollaries 

A number of results follow right away from the existence of the embedding map (B exhibited 
in the previous section. We shall describe those next. 

(1) In summary, if C'j' is an eigenvector of weight w with Laplace-eigenvalue s(s — 1) 
for s = m — 2m and m = 1, 2 ,..., [(m — l)/2 ], then for p > 1, the function is an 
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eigenvector of weight w + 2^ with the same eigenvalue s(s — 1). It is manifest by inspection 
that the map C is injective. Thus, we have 


V(w + 2,s) = e(V(w,s)) 

dim V{w + 2, s) = dim V(tc, s) 


s = w — 2m 
m = 1,2,... 



(3.49) 


(2) Given the first emergence of the eigenvalue s(s — 1) at weight s = w — 2 hj the 
eigenvector combinations Cj, the eigenvectors corresponding to the same eigenvalue s(s —1) 
and higher weight w = s — 2 + 2fi are obtained by fE^{Cj). To prove the structure of the 
spectrum, it will thus suffice to prove that. 


dim V(tc, w — 2) 



(3.50) 


MAPLE calculations up to weight w = 34 conhrm this assertion holds. A complete proof of 
this structure will be given in section 4. 

(3) The map (£ may be expressed in terms of the mutually commuting shift operators 
Ta,Ti,, and acting on the variables a, b, c and on the functions Ca,b,c by, 

TaCa,b,c = Ca+l,b,c [Ta, d] = Ta (3.51) 

with analogous relations for Tb and Tc, so that we have 


(B = obTaTb + bcThTc + caTcTa 


(3.52) 


The fi-th power of the map € is obtained via the trinomial formula and we hnd. 


«'‘(Ca,v)(T) = E 


fi\T{a + fd- Ai 2 )r(& + fdi+ Ai 2 )r(c + fd- fii) 
„ „ To r(/i - /Xi - /i 2 + 1 ) r(a) T{b) r(c) 

X Ca+^—ii2,b+fii+fi2,c+fj,—^ii {t) 


(3.53) 


The predictions of this formula have been tested against explicit MAPLE calculations at 
weights up to tc = 34, and agree. We shall limit ourselves here to providing examples up to 
weight 11, and make contact with the examples given to this weight order in section 3.5. 


3.8.1 Illustrations up to eigenvalue 12 

At eigenvalue 0, with multiplicity 1 for all odd weights, we hnd, 

^(^* 1 , 1 , 1 ) = 3(72,2,1 (3.54) 

^(C*2,2,l) = 4(73,3,1 + 4(73,2,2 

^(^*3,3,1 + ^*3,2,2) = 9(74,4,1 + 18(74,3,2 + 4(73,3,3 
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Comparison of the right hand sides with the entries of (3.33) shows agreement on the top 
three lines, while the last line above predicts the correct linear combination for weight tc = 13. 
Similar comparisons at higher eigenvalue further conhrm the validity of the general formula. 
At eigenvalue 2, we have, 


^(C*2,l,l) = 4(73^2,1 + ^*2,2,2 (3.55) 

^(46*3,2,1 + ^*2,2,2) = 4(6(74,34 + 3^*4,2,2 + 5(73,3,2) 

^( 6 ( 74 , 3,1 + 3 ( 74 , 2,2 + 5 ( 73 , 3 , 2 ) = 9 ( 8 ( 75 , 4,1 + 8 ( 75 , 3,2 + 7 ( 74 , 2,2 + 8 ( 74 , 3 , 3 ) 

which agrees with (3.35). Skipping eigenvalue 6, which is analogous to eigenvalue 2, we have 
for eigenvalue 12 where the hrst degeneracy occurs. 

For eigenvalue 12, we shall just list here the map between weight 6 and weight 8, 
^(6(73,2,1 — (72,2,2) = 18(2(74,3,1 + (74,2,2) 

^(6(74,1,1 + ( 72 , 2 , 2 ) = 6(8(75,2,1 + (74,2,2 + 2 ( 73 , 3 , 2 ) (3.56) 


which again agrees with (3.37). 


3.9 Explicit algebraic relations for s = 1 and for all odd weights 

It was established earlier that the zero-eigenvalue equation on the functions Ca,b,c for any 
odd weight w = a + b + a leads to an algebraic equation between the (7-functions and the 
Eisenstein series E^. The simplest one of these relations was derived in (3.11), and a list of 
the hrst cases up to weight 10 was given in (3.33). We shall now present an explicit formula 
valid for all weights. 

Parametrizing the weight by tc = 2/i + 3 for /i a positive integer, we may obtain the 
linear combination of (7-functions at weight w which give the eigenfunction of the Laplacian 
for weight 0 by applying the embedding operator to the function with zero eigenvalue 
at weight 3, namely (7i,i,i with the help of the general formula (3.53). The right hand side 
of the zero eigenvalue equation at weight w is proportional to E^. Therefore, the Laplace- 
eigenvalue equation for zero eigenvalue can be integrated to an algebraic equation at any 
odd weight w, and this equation is given by. 


E = fMr) + (3.57) 


/7l ,7*2=0 


/ii! /i2! (+ — +1 — +2)! 


The summand vanishes for /xi + /i 2 > + so that the sum truncates, as expected. Here, 
and Qw are constants which depend only on w, and may be determined from the asymptotic 
behavior near the cusp; their values for weight up to 9 are given in (3.34). 
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4 The generating function for Ca,i,c 


It is natural to introduce a generating function for the modular functions Ca,b,c with a, 6 , c > 1 
in terms of three independent real variables 

OO 

>V(tl,t2,t3|'r) = CafiA'T) (4-1) 

a,6,c=l 


The lowest order term is independent of and reduces to Ci^iAA- Since the Ca,b,cA) 

are S'L(2, Z)-invariant modular functions of r, so will W be provided the variables 
transform trivially under 5'L(2,Z), which is what we shall henceforth assume. Using the 
explicit expression for Ca^AA given in terms of its dehning sum in (3.1), we may carry out 
the corresponding geometric series, and we hud the following expression for W, 


>V(ti,f2U3|T) 


(mr,nr)^(0,0) 


^n,0 

r=l 


T2 


nlnirT + nA — trT2 


(4.2) 


Recall that we use the notation m = mi+m 2 + ms and n = ni+n 2 + ns for the components 
of total momentum, and that the Kronecker S factors enforce the vanishing of m and n. 

Setting —tr = would give the 2-loop sunset diagram on a torus with Euclidean 
signature and modulus r for three free scalar Green functions with masses Mf, M|, and M|. 
For fixed r = Ti + iT 2 with r 2 > 0, the above sum is absolutely convergent for small enough 
ti,t 2 ,ts. This may be established as follows. Contributions to the sum with 7 ^ 0 are 
absolutely convergent provided \tr\ < 71 T 2 , while those with mr = 0 are absolutely convergent 
provided |G| < 7 r/r 2 . Thus we have convergence when |U| is less than both nT 2 and 7 r/r 2 . 
This convergence is uniform in Ti for hxed T 2 , but fails to be uniform in r throughout A^i, 
as the sum diverges for T 2 ^ 00 . 


4.1 Differential equation for the generating function 

We shall now translate the algebraic representation of the Laplace-Beltrami operator A on 
the functions Ca,bAA given in (3.6) into a representation on the function W by a linear 
differential operator in ti,t 2 ,ts. To do so, we use the algebraic representation of A on Cafi,c 
to compute its action on W. This may be done term by term in (3.6). A Erst example is 
given by the first term on the left side of (3.6), for which we have, 

CxD 

a{a - 1 ) Ca,b,c = hdAhyV) (4.3) 

a,6,c=l 

where we use the notation dr = d/dtr for r = 1,2,3. Clearly, the terms given by b{b — 1) 
and c(c — 1) are obtained by permutations of ti, t 2 , ts- A second example is provided by the 
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first term on the right side of (3.6), for which we have, 


l,6+l,c EE tli-HrHa + l){b-l)Ca,b,c (4.4) 

a,6,c=l a=0 6,c=l 

This expression is compared with the action of the following differential operator on W, 

CxD 

d,d2{tlW)= Y, (a + l)(&-l)t?4"2tr'C'aAc (4.5) 

a,b,c=l 

The two expressions match, except for the contribution from a = 0, which may be expressed 
purely in terms of Eisenstein series with the help of the first equation in (3.8). Note that it 
will involve the regularized quantity E^, which in the end will consistently cancel out of all 
results, as we shall show explicitly below. The second term on the first line of (3.6) as well 
as the first two terms on the second and third lines may all be obtained by permutations of 
ti,t 2 , ts- The same analysis applies to the remaining terms on the right side of (3.6). Putting 
all together, we have. 


A W - £2 W = 

where we have introduced the following operators. 


(4.6) 


D = tidi+ t2d2 + hds 

= (2)^ + 2}) + [tl + ^2 + ^3 — ‘2tit2 — 2t2^3 ~ (^did2 + did^ + d2d^ (4.7) 
The right hand side function TZ is given by. 


n = 


oo 

E 

a,b=0 


oo 

(ti ^2 +^2^3+ ^3 ti j + Yj (^1 ^2 ^3 + ^2 ^3 + ^3 ^1 ^2 j (4.8) 


a,6=0 


in terms of the following combinations of Eisenstein series, 

TT-af, = 3a(& + l)Ea+iEb ^2 + 3(a + l)bEa+ 2 Eb+i + {2 — a — b — 4ab)Ea+b+3 
'^ab ~ (ib(^Ea+2Eb+2 ~ Ea+b+4^ (4-9) 


Note that and 1Z\^ are symmetric in a, b, and that they never involve E^ in view of the 
fact that the corresponding coefficients that enter 77°^ for a = 0 or 6 = 0 cancel. 
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4.2 Adapted coordinates and realization of ©3 

The operators D and 2“^ are manifestly invariant under rotations of the variables ti,t 2 ,t 3 
which leave the direction (1,1,1) invariant. Cyclic permutations of these variables then 
correspond to 120 degree rotations around this axis, while permutations with odd signature 
correspond to reflections through any plane that contains this axis. We shall now make a 
linear change of variables in which these transformations are realized naturally, 


U — (ti + t 2 + ts) 
V3 


V = 


V = 


(ti + et2 + 

{tl + £^t2 + ets) 




In terms of these variables, we have, 

^2 — 2tit2 — 2t2t3 — 


8182 + 8183 + 8283 = 


—u^ + 2vv 
8l — 28„8r, 


and the operators become, 


( 4 . 10 ) 


( 4 . 11 ) 


2 ) = u8u + v8v + v8y 

2 ^ = £+£_ + £_£+-£2 ( 4 . 12 ) 

where the first order operators £±,£3 are expressed as follows. 


£+ = u8y + v8u 

[£0, £+] — +f£+ 


£_ = u8y + v8u 

[£o,£_] = -f£_ 


£0 = iv8y — W8y 

[£+, £_] = -i£o 

(4.13) 

The algebra is that of S0{2,1) = S'0(2,1;M), and £^ is its quadratic Casimir operator. 
Both combinations in (4.11) are invariant under SO{2, 1). 

Under the permutation group ©3, the variables 

u, V and v transform as. 


{123){u,v,v) = {u,v,v) 

(132)(m, n, h) = {u,v,v) 


(231)(M,n,h) = {u,e‘^v,ev) 

(213)(M,n,h) = {u,ev,e‘^v) 


(312)(M,n,h) = {u,ev,e‘^v) 

(321)(M,n,h) = {u,e‘^v,ev) 

(4.14) 


could stay with real coordinates, and decompose v into real coordinates x,y hy v = {x + iy)j'/^. We 
would then have £ = udu + xdx + ydy as well as 2^ = £g with £i = udx + xdu, £2 = udy + ydu 

and £0 = xdy — ydx, so that each one of these operators is manifestly real, thereby justifying reference to 
the algebra SO{2, 1;M). The operators £1 and £2 are related to £± by £± = (£1 ± 122)1'^- 
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The generating function W is real; it is a sum of monomials of the form u°'v^v'^; and it is 
invariant under the group © 3 . We shall now use this information to constrain the form of 
W when expressed in terms of the variables u, v, v. To do so, we associate to an arbitrary 
monomial u°‘v^v'^ of weight w = q;+/ 9 + 7+3 a unique corresponding © 3 -invariant polynomial 
obtained by summing over the six images of the monomial under permutations © 3 , 

©3 : (v^v^ + (l + + e2(/3-7) j (4.I5) 

Clearly, the last parenthesis vanishes when — 7 ^ 0 (mod 3), in which case the monomial 
produces no © 3 -invariant contribution to W. Therefore, all the non-trivial © 3 -invariant 
contributions of weight tc to W must take the form, 

^ ^ 3 p j w = s + 2m (4.16) 

where m > 1 and p > 0 with s — 1 > 3p. Changing basis we can equivalently work with the 
following set of functions. 


_ 2^;^)’^-! (^^;3p ^ -3p^ 

This basis has the advantage that the combination u'^ — 2vv is invariant under SO{2, 1), so 
that — 2vv) = £o('W^ — 2vv) = 0 . 

The image of a monomial under £0 is not invariant under 

© 3 , so that £0 does not have a good action on the © 3 -invariant functions W. However, the 
action of the operator £g is well-dehned on this space. Thus, we shall choose to simultane¬ 
ously diagonalize £>, £^, £q. 

4.3 Diagonalizing X), i^o 

We choose to simultaneously diagonalize the three commuting operators, T>, £^, and £q, and 
denote the corresponding eigenvalues respectively by {w — 3), s{s — 1), and —for p > 0, 

2) >Vu,;s;p = (m - 3) >Vu,;s;p 
-2^ >V^;s;p = s(s - 1) >V^;^;p 

= -9p^W^.,s;p (4.18) 

The first and last eigenvalue equations in (4.18) are manifestly satished by construction. 
The representation theory of SO{2, 1) will relate these eigenvalues, including s, as we shall 
establish shortly. Weight w signihes that the eigenfunction Ww,s-,p is a homogeneous polyno¬ 
mial of degree m — 3 in the coordinates (ti, ^ 2 , ts), or equivalently in {u, v, v). Therefore, the 
corresponding representation of *S'0(2,1) will be hnite-dimensional. 
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The eigenfunctions VVw-s-,p take the form, 


m+ 

fc=m- 






(4.19) 


for as yet undetermined real coefficients J^w,s-,p{k). The integers m± are subject to the con¬ 
ditions 1 < m_ < m+ < tc — 3p — 1. We shall assume that rri- is the lowest value for k with 
non-zero contribution, while m+ is the highest value, so that J^«,,s,p(m±) 7 ^ 0 , but otherwise 
m± remain to be determined. It is straightforward to compute the operator £? on these 
functions. The eigenvalue equation for £? in (4.18) then reduces to the following equations. 
For m_ < k < m+ we have. 


{w — s — 2k){s + w — 2k — l)J^w,s,p{k) 

= {w-2k-3p + 1)(m; -2k- 3p)J'^,,,p(A; - 1) (4.20) 


while for the extremities, we have the conditions (using the assumption J^u,,s,p(tn±) 7 ^ 0 ), 

(tc — s — 2m._)(s -I- tc — 2m_ — 1) = 0 
(w — 2m+— 3p — l)(tc — 2m.+— 3p — 2) = 0 (4-21) 

The hrst equation forces two solutions for s which are equivalent as they yield the same 
eigenvalue s(s — 1) for We shall make the choice s = w — 2m._. The solution of the 
second equation is easily identihed as well, and assembling the full solution, we have. 


s = w — 2m_ 


m+ 


tc — 3p — 1 
2 


(4.22) 


with m_ > 1 , s > 1 , and p > 0 . 


4.4 Proof of Theorem 1 


The structure of the spectrum of the Laplace operator, which had been inferred on the basis 
of MAPLE calculations in Theorem 1 of section 3.5, may now be proven with the help of 
the above results on the structure of the generating function W. In fact, the structure may 
be read off from the relations of (4.22). 

For given weight w, the allowed eigenvalues s(s — 1) of the Laplacian are governed by 
s = tc — 2m_ with m_ = 1, • • •, [(tc — l)/2]. Identifying m_ with m in Theorem 1, we see that 
we have the same set of eigenvalues for given weight w. Their degeneracy may be read off 
from (4.22) as well. Combining the condition m_ < m+ with the second equation in (4.22), 
and expressing w in terms of s and m_ gives. 


m_ < 


s -I- 2m._ — 3p — 1 
2 


(4.23) 
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Subtracting m_ on both sides then gives a condition which equivalent to 0 < 3p < s — 1, 
which in turn implies that the degeneracy of the eigenvalue s(s — 1) is given by [(s + 2)/3], 
as indeed announced in Theorem 1. 

In addition, we obtain an explicit form for the eigenfunctions given by their coefficients, 




T{k+ T{k+ 

T{k + T {k- "" 


(4.24) 


where is an arbitrary constant, and k is in the range m_ < k < m+. 

Note that the map (£ of section 3.7 corresponds to multiplication by a factor {v? — 2vv), 
which leaves the S'0(2,1) representation unchanged, but increases the weight by 2. 


4.5 Illustration for low weight w 


Taking the normalization J7«;s;p(iTi-^ 
W’^;s;p(n, n, h) up to weight w = 6, 

) = 1, we may now write 

down the basis functions 

y^3;l;0 

= 1 


y^4;2;0 

= u 


kV5;3;0 

= — — 2vv) 


>V5;l;0 

= — 2vv 


>V6;4;l 

= + v^ 


kV6;4;0 

3 

= - uiu^ — 2vv) 

5 


kV6;2;0 

= u{v? — 2vv) 

(4.25) 


The expansion coefficients of W in this basis are the eigenfunctions referred to in 

(1.5) of the Introdnction, so that the generating fnnction may be written as follows. 


>V(ti,t2,t3|'r) = 2m;p(^5 1 ^w;w — 2m;p (r) (4.26) 

ui=3 m=l p=0 

where the eigenvalue s is given by the combination s = w — 2m. The functions £u,;s,p(t) may 
be read off by decomposing W onto the basis functions Ww-s-p, and we hnd. 


^3;l;0 

^4;2;0 

^5;3;0 


— C*!,!,! 


= \/3C*2,1,1 


= 4(2C'3,2,1 

+ Q 
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^ 5 ; 1;0 

^ 6 ; 4;1 

^ 6 ; 4;0 

^ 6 ; 2;0 


2 <^ 2 , 2,1 

^ 3 ( 744,1 — 3 ( 73 , 2,1 + ^* 2 , 2 , 2 ^ 
^ 30 ( 74 , 1,1 + 24 ( 73 , 2,1 + < 72 , 2 , 2 ^ 

10 


4 ( 73 , 2,1 + < 72 , 2 , 2 ^ 


(4.27) 


Comparing the entries with s = 1,2,3 with the expression for s = 1 given in (3.33), s = 2 
given in (3.35), and s = 3 given in (3.36), we hnd perfect agreement. For s = 4 we hnd perfect 
agreement with (3.37) upon taking linear combinations which had been left unspecihed in 
(3.37), but which are here dictated by the fact that the eigenfunctions also diagonalize £q. 


5 Modular functions in weight-4 interactions 

In this section we shall consider the weight-4 Feynman diagrams that contribute to j42’°^(t). 
Their integrals over A4i will determines the coefficient of the one-loop interaction. 

The weight-4 diagrams that will enter this calculation are 71i,i,i,i, and = (7i,i,i,i. 
The precise rational coefficients with which they contribute to j 7 ,o) listed below in 

(7.12). The modular function D 2 equals the Eisenstein series E 2 , while Di i i i equals but 
the modular function ZI 4 = (7i, 1 , 1,1 has a novel form beyond the framework of the functions 
Ca, 6 ,c studied in the preceding section, though it does £t into the framework of the functions 
defined in (1.1) of the Introduction for the case for p = 4. 

We have not found a useful closed algebraic representation for the Laplace operator acting 
on the modular functions (7ai,...,ap(T) for p > 4. Thus, we do not have generalization of the 
powerful tool (3.6) at our disposal. It would be very valuable, of course, to find such an 
algebraic representation, if it exists at all. Instead, our procedure will be based on expanding 
the multiple sums that define the modular function in the vicinity of the cusp in the large 
T 2 limit, and using the resulting expansion to conjecture the form of the Laplace operator 
acting on D/^. In section 6 we shall extend this procedure to modular functions of weight 5. 

5.1 Asymptotic expansion near the cusp r —ioo 

We shall begin by producing the general form of the expansion near the cusp r —)■ ioo. 
To minimize the occurrence of factors of tt throughout, it will be convenient to express the 
asymptotic expansion of the modular functions (7, D, and E in terms of the parameter, 

y = 7 rr 2 (5.1) 
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The asymptotic expansion near the cnsp y —)■ cx) of a modnlar fnnction D (and similarly of 
the fnnctions C and E) takes on the following form, 


D{t) = (5.2) 

k,k=0 

where q = and each coefficient fnnction T>^^’^\y) consists of a hnite snm of integer 
powers of y, namely it is a Lanrent polynomial in y of hnite degree. The series is absolntely 
convergent at any point in the interior of A^i. Our conjectured relations will be verihed by 
comparing the pure power-behaved terms for /c = A; = 0 as well as the leading exponential 
corrections for A; -|- A; = 1. For this purpose, it will be convenient to have also the expansion 
of the Laplacian of D{t) available, which is given by, 

CXD 

AD(r) = '^{AVY’^’^\y)q^q^ 

k,k=0 

{AVf^’^\y) = y^dlV^^'^\y)-A{k + k)y^dyV^^’^\y) + lQkky^V^^'^\y) (5.3) 

The Fourier transform of the function D{t) in the variable Ti at hxed value of y is then 

obtained by the following inhnite sums, 

« OO 

/ dTi D(r) = I,(5,4) 

k,k=l 

An important special case is the expansion is for the Eisenstein series En itself given in 
(3.24) where the iV-th Fourier modes gets a contribution from the single power q^. More 
generally an inhnite number of coefficient functions (with k — k = N) contribute to a 

given mode. These appear to correspond to the ehect of an inhnite number of world-sheet 
“instanton-anti-instanton” pairs (an interpretation which deserves further study). 

5.2 The case 02,1,1 = ^ 2 , 1,1 

It will prove of interest to hrst consider a known example and examine the properties of 
the weight-4 modular function O,!,! = £> 2 , 1,1 with the help of the asymptotic expansion 
procedure discussed above. Its Feynman diagram is given as follows. 
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The function 6 * 2 ,i,i satisfies the Laplace equation (3.15), but will not arise directly in the 
expression for j( 2 ,o) given in (7.12). We can usefully rewrite the associated Laplace- 

eigenvalue equation as, 


(A-2) 




(5.5) 


The power-behaved terms in the expansion of 6 * 2 ,i,i near the cusp y = oo can be 

obtained from this equation in a straightforward manner by expanding the source term and 
matching the power-behaved terms. The y~^ and y"^ terms are in the kernel of (A — 2 ) 
and are determined by enforcing appropriate boundary conditions, or alternatively may be 
determined by carrying out explicitly the multiple sum which dehnes ( 7 * 2 ,i,i- One hnds, 

(0,0) _ ^ ^ ^ C) 

14175 ^ 45 ^ 12y 4y^ ^ 16y^ ^ ’ 

The coefficient of the leading exponential correction terms may be obtained by con¬ 

sidering the first Fourier mode of (5.5) truncated to its lowest order exponential part, and 
doing so results in the expression. 


( 1 , 0 ) _ y 1 11 C(3) 9 9 

45 3 12y 2y^ 8y^ 16y^ 


(5.7) 


Here we have already enforced the condition of power growth near the cusp to set to zero a 
contribution which would grow exponentially at the cusp. A term proportional to 1 -|- l/{2y) 
in the formula above is in the kernel of (A — 2 ) and may be determined by boundary 
conditions. Alternatively, this expansion may be obtained by direct approximation of the 
multiple sums that dehne ( 7 * 2 ,i,i- The power-behaved terms were analyzed in this manner 
in [25], and the two approaches are found to agree. Here we have also determined the leading 
exponential correction. 

The above exercise is somewhat redundant in the case of ( 7 * 2 ,i,i since we have derived the 
Laplace equation that it satishes. However, the direct expansion of other C- and H-functions 
will prove signihcant in the sequel where the Laplace equations are conjectural. We now turn 
to consider the function which will enter in the expression for j( 2 ,o) (7.12). 


5.3 The case TI 4 = Ci, 1 , 1,1 

The Feynman diagram representation of II 4 is as follows, 

•ATEA* = D,(.Ti 
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Although we have not determined the functional form of D^lr) directly from its expression 
as a multiple sum, we will arrive at a convincing conjecture for its form, based on examining 
its behavior near the large y cusp. The power-behaved part of the expansion was given 
already in equation (B.9) of [25], 

,5.8) 

945 6 y 

The leading exponential corrections may be evaluated in a similar manner, as will be detailed 
in appendix A, and we find. 




(») 


4i/V 2y 4 12C(3) 

15 3 V y 


6C(3) 99 

y2 2|/2 4y3 


(5.9) 


Since this function is manifestly real, we have vf’^\y) = V^^’^\y). It may be verified by 
explicit calculation that the above pure power part and leading exponential corrections in 
the expansion of 1^4 satisfy the pure power part and leading exponential corrections of the 
following equation. 


(A-2) 



-2AE^ 


(5.10) 


Based on this evidence, along with its overall modularity and weight, we conjecture that this 
equation holds true exactly for all r in the upper half plane. 

Assuming the conjectured equation (5.10) holds true exactly leads to a new algebraic 
equation involving D^. Indeed, eliminating the term between the right hand side of 
(5.10) and the right hand side of (5.5) leads to the following Laplace-eigenvalue equation, 

(A - 2) {Di - 24^2,1,1 + 18^4 - 3E^) = 0 (5.11) 

Given that the argument of the above equation has power behavior near the cusp, we know 
that its solution must be proportional to E 2 . Matching the pure power behavior suffices, 
however, to show that the E 2 contribution is in fact absent, so that the unique solution is 
given by the algebraic relation, 

D 4 = 2402,1,1 - I 8 E 4 + 3Ei (5.12) 

In the next section we shall produce analogous conjectures for all the weight-5 Feynman 
diagrams by expressing them in terms of modular functions Ca,b,c of weight a -|- 5 -|- c = 5 
and Eisenstein series. 
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6 Modular functions in weight-5 interactions 

We will now consider the weight -5 modular functions 1,1,1, -Da,!,! = (^2,1,1,1, D2,2,i 

and D2,1,1,1 = Da,1,1 that contribute non-trivially to and will contribute to the interac¬ 
tion of order As in the last section, we will not determine these functions directly 

from their dehnition in terms of integrated powers of Green functions, but we will study 
their asymptotic behaviour near the cusp. In this manner we will be led to relationships 
that will express the functions D5, D3,i,i and D2,2,i in terms of (^3,1,1, which satishes the 
Laplace equation ( 3 . 19 ). We will now consider each of these in turn. 

6.1 The case -D2,1,1,1 = C3,i,i 

The Feynman diagram representation for 6*3,1,1 is as follows. 



= a 


3,1,1 


:r) 


This modular function satisfies the Laplace-eigenvalue equation ( 3 . 19 ), which may be re¬ 
written in the following form. 


(A-6) 



43 


D5 -|- 


60 


— 4 D 2 D 3 


( 6 . 1 ) 


For later reference we give the pure power part (which was also given in (B. 53 ) of [ 25 ]), 

2y^ 


cSb) = 


155925 


2C(3)y^ _ C(5) 707) _ C(3)C(5) 43C(9) 

' ^ A -I r\ ' -I O ' 


945 


180 16 i /2 


2i/3 


64i/^ 


( 6 . 2 ) 


Furthermore, the leading exponential correction term can be extracted directly from the 
dehnition of 6*3,1,! as a multiple sum, or by considering the hrst mode of the equation (6.1). 
The results agree and take the form. 


^(1,0). 5 _ ^ 2m 

^3,1 Ay) - 945 + 45+ 3 


45 y 


43 C(3) C(5) 
24i/ 2i/3 2 i/3 


43 


43 


32i/3 


(36 3 

+ C5 4 H-1 —- 

y y 


( 6 . 3 ) 


The constant C5 is not determined by the mode expansion of ( 6 . 1 ), but may be evaluated 
directly from the sum representation of 6*3,1,1 and we hud. 


1 19 


(6.4) 


Since 63^+ is real, we have 63[^(]^/(i/) = 63[f[’/(i/). 


( 0 , 1 ), 


Fi,o), 
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6.2 The case = 01,1,1,1,1 

The Feynman diagram representation of D 5 is given as follows, 



Here again we have not obtained a Laplace equation for directly from its definition as 
a multiple sum. However, as with ^ 4 , we conjecture a Laplace equation for based on 
its asymptotic expansion near the cusp. The power part of the expansion in this limit was 
given in (B.IO) of [25], and takes the form. 


pr\») = 


V ^ 10y"C(3) ^ 95C(5) ^ 10C(3)^ 


18711 
^ 105C(7) 


27 Q y 

45C(3)C(5) , 225C(9) 


4|/2 2y^ 

while the leading exponential correction takes the form. 


16?/"* 




255 


8 y^ lOy^ 5y 535 80C 3 

189 27 3 ^ 6 y 2y 


15C(3) 15C(5) 

o o 


y 




45 

y2 


45C(3) 45C(5) 225 225 


2?/3 


2?/3 


8?/3 l&y^ 


Q 

+60 C 5 ( 4 + - + — 

y y 


(6.5) 


( 6 . 6 ) 


where the constant C 5 was determined from the sum representation of in appendix A, and 
is found to be given by (6.4). Since V^’^\y) is manifestly real, we have V^^’^\y) = V^’^\y). 

If we now assume that is related to a superposition of the other weight-5 functions 
that arise in the four-graviton amplitude, then one is led to conjecture the following Laplace- 
eigenvalue relation. 


(A - 6 ) {D 5 - lOEsCi,!,!) = 360^5 - 24 OE 2 E 3 - 90C(5) (6.7) 

Eliminating the product E^E^ between this equation and the Laplace-eigenvalue equation 
for 6 * 3 , 1,1 established in ( 6 . 1 ), and moving E^ under the operator A — 6 using the eigenvalue 
equation (A — 20 ) 7^5 = 0, we obtain the following relation. 


(A - 6 ) - hoc's, 1,1 - IOE 3 E 2 + 48^5 - 16C(5)) = 0 ( 6 . 8 ) 

One verihes that within the pure power and leading exponential correction approximation, 
the above equation is indeed satisfied. With this conjecture (and noting that the large y 
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expansion fixes the coefficient of the solntion of the homogeneous equation to be zero) we 
conclude that is given in terms of -D 2 , 1 , 1,1 = ^* 3 , 1,1 by the following algebraic relation, 

Ds = 60^3, 1,1 + 10^2^1,1,1 - 48^5 + 16C(5) (6.9) 

which bears a striking similarity with the corresponding equation for in (5.12). 

6.3 The case = C 2 , 1 , 1,1 



Once again, we have evaluated the complete pure power part for this function in the 
large-y limit (which was partially given in (B.34) of [25]). This has the form, 

^ Jiq y\(3) 11C(5) 105C(7) _ 3C(3)C(5) 81C(9) 

3 . 1,1 W; 22275 ^ 45 60 2tf< My* 

The asymptotic expansion leads to strong constraints on the possible form of the modular 
function, and we are led to conjecture the following algebraic relation involving £* 3 , 1 ,i, 

40T)3,i,i = 30003 , 1,1 + I 2 OO 2 O 3 - 276^5 + 7C(5) (6.11) 

We have verified this equation only to pure power order. 

6.4 The case /^ 2 , 2 ,i 

The Feynman diagram representation od 7^2,2 ,1 is given as follows, 

= -02,2,1 ('r) 




This is an example of a modular function that is not of the form C'a 2 ,...,ap for any value of p. 
However, we will see that its asymptotic behavior leads us to conjecture a relation again to 
the modular function 02 , 1 , 1,1 = 03 , 1,1 of weight 5, along with Eisenstein series. In this case 
the expansion in the large y limit gives (noting some errors in the terms in (B.35) of [25]), 


8i/5 4C(3 )i/^ 13C(5) 7C(7) C(3)C(5) 9C(9) 

467775 945 45 8i/2 Sy^ 


( 6 . 12 ) 
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The expansion leads to conjecture the following Laplace-eigenvalue equation, 


(A-6) 





—8E^E2 


(6.13) 


Combining (6.1) with (6.13) leads to the following conjectured algebraic relation 


10/12,2,1 — 20(7344 ~ 4/?5 + 3C(5) 


(6.14) 


This equation may be used to predict the leading exponential corrections (y) = iv )) 
and we hnd, 


2 , 2,1 y y2 y2 y3 yZ ^yZ ^yi 


(6.15) 


By directly evaluating this sum, we have checked in appendix B that all terms indeed match. 

This concludes our study of the weight-5 modular functions that contribute to the higher 
derivative interaction and we shall now proceed in the next section to providing 

further numerical evidence for the conjectured relations, and to evaluating the integrals over 
All of all these contributions in section 7. 


7 Low order coefficients in the genus-one expansion 

We will now obtain the integrated coefficients of some low order terms in the low-energy 
expansion of the amplitude (2.3). The coefficient of a term of order af crl is obtained by 
integrating dehned in (2.11), over the fundamental domain, Ali. Since j^^’'^^(t) is 

C>(r|) as T 2 —)■ cxD, with a > 1, such integrals diverge. Our discussion will begin with a brief 
review of treatment of these divergences, which are related to the occurrence of non-analytic 
threshold terms in the expansion. We will also review the explicit determination of the 
coefficients of terms with weights w < 3 before determining the coefficients of the weight 4 
and weight 5 interactions^^. 

7.1 Separation of analytic and non-analytic terms 

An important new issue appears in considering the low energy expansion of the amplitude 
beyond the tree-level approximation. This concerns the appearance of physical threshold 
singularities that are branch cuts in the Mandelstam variables, and correspond to non-local 


weight-4 and weight-5 results correct the discussion in the earlier literature [4, 25] where only 
the asymptotic behavior at large y of the relevant modular functions was evaluated and the values of the 
coefficients that follow from integration over r were consequently not correctly determined. 
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contributions to the effective action originating from integrating out massless modes. In 
order to make sense of the power series expansion (or local terms in the effective action) 
it is necessary to separate the analytic and non-analytic parts of the amplitude, which is 
straightforward up to the orders we are considering here. This is achieved by splitting the 
fundamental integration region in (2.1) into two sub-domains [25], 

M, = M,lIJM,r (7.1) 


where the first component is the cutoff fundamental domain. 


Mil = {r ; 



t\ >1, r 2 < L} 


(7.2) 


and L is a positive real parameter which may be taken to be large L 3> 1. The remainder 
Mir of the fundamental domain A^i is the rectangular region that bounds the cusp, namely 
Mir = {t ; |ti| < T 2 > L}. Accordingly, the one-loop amplitude Jii\ei,ki) is given by 
the sum of two parts. 


k,) = k,) + .45“'"“-“"(e,, L) 

The analytic part of the amplitude, .45“'“"(e,, /t,), given by, 


(7,3) 


^( 4 ) an 


(ei,/Ci) 


27r k-Iq 77^ 



d/ii Bi{s, t, u\t) 


(7.4) 


admits an expansion in powers of <72 and < 73 , in view of the fact that the integrand Bi{s, t, u\t) 
admits a series expansion in <72 and < 73 , 


Bi{s,t,u\T) = al (7.5) 

p,q=0 


which is uniform in r throughout the region Mil- In particular, the integrals of the coefficient 
functions over Mil are convergent, but do depend upon the regulator L. 

The non-analytic part of the amplitude ki), is given by. 


^{ 4 )non = 271 f d^i Bi{s,t, u\t) ( 7 . 6 ) 

•J Mm 

Although the expansion of (7.5) continues to hold pointwise for r in the interior of Adi, it 
does not hold uniformly in Mir, so that the integral of Bi over Mir will fail to be analytic 
in s,t,u near the origin. In particular, the integrals of the coefficient functions j*'^’'^^(t) 
over Mir will diverge for large enough p, g, further signaling the breakdown of analyticity 
of Thus, in the region Mir, the expansion (7.5) is invalid; the integral 
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needs to be treated non-perturbatively in s, t, u. The resulting value of “"(ej, ki) will 

depend upon L, but the dependence on L in the full amplitude will, of course, cancel. 

Concretely, the low-energy behavior of ki) is determined by expanding the 

integrand around T 2 = 00 , thereby generating threshold singularities in s, t, u associated with 
the propagation of massless intermediate states. The hrst of these is the supergravity thresh¬ 
old of the symbolic form Sij log(sjj), which can also be determined directly from dimensional 
regularization of maximal supergravity in 10 — e dimensions. However, in the string theory 
context there are further “stringy” thresholds, which were discussed in detail in [25]. The 
hrst of these arises at order log(sij). 

The analytic part of the low energy expansion of the one-loop four-graviton amplitude is 
obtained by integrating (7.5) term by term over the cutoff fundamental domain, 

(7.7) 

JMil 

In the limit L 3> 1 this gives an expression of the form, 

(L) = + c log(L A) + + 0{l/L) (7.8) 

i 

where a^, c, A and r* are constants, with Ti > 0. This is a sum of positive powers of L, 
together with a log L and a constant term. 

The L-dependent terms are cancelled by corresponding terms from the non-analytic 
part, leaving the sum of clog A and The presence of clog A signals the presence 

of a threshold contribution to the non-analytic part of the amplitude of the symbolic form 

log(sjj//i) while is the coefficient of the interaction of order af Cg A re- 

dehnition of the scale, /i —)■ e~^ fi of the logarithm leads to a redehnition of the analytic term 
by the addition of pa^ Cg 77^. This generally implies an ambiguity in the dehnition of certain 
analytic terms, as we will see in the weight-4 example later in this section. A more detailed 
discussion with explicit examples is given in [25]. The fact that integrals over a fundamental 
domain of functions of moderate growth, such as j*'^’'^^(t), can be interpreted as sensible 
hnite quantities has a respectable mathematical pedigree, as discussed in [43]. An example 
of immediate relevance is that with this interpretation the integral of a non-holomorphic 
Eisenstein series, which vanishes, dpiEg = 0, if s > 1, since = 2sC(2 s)L* 1 , 

which has no L-independent piece. The resulting expansion of the amplitude takes the form 
given in (2.12). 

The modular invariant coefficient functions jAp-?) dehned by (2.11) are given by sums of 
Feynman diagrams with coefficients that are simple combinatorial factors that arise from 
expanding the right-hand side of (2.10). The resultant sums of diagrams that arise at low 
orders in the expansion were presented in appendix C of [25], and are summarized here for 


50 



reference in the following sections: 



1 

(7.9) 

4(1.0) = 

D 2 

(7.10) 

4(0.1) ^ 

^(8-Di,i,i -I- 2D3) 

(7.11) 

4(2.0) ^ 

^(^4 -|- 9D2 -|- 6D1 1 1 1) 

(7.12) 

4(1.1) ^ 

— (2D5 -|- 32D3 11 — 24D2 2,1 + 24D2 111 

0! 

— 48 Di 11 i-i -|- 28D3D2 -|- 96 Di 11D2) 

(7.13) 


All the D functions that enter into these expressions up to weight 5 have been determined 
in the previous two sections, so we will now evaluate the integral of for all w = 

2p + 3g < 5. 

7.2 Interaction coefficients with weights < 3 

The integration of the modular invariant coefficient functions up to weight 3 over a funda¬ 
mental domain were determined in [4], The results are summarized as follows. 

The weight-0 diagram has no propagators and Do = 1, so the coefficient of is given 
by the volume of the fundamental r domain. In other words, 

5(0,0) ^ I (714) 

The only weight-1 diagram is the zero mode of a single propagator, which corresponds to 
the vanishing of D^TZ"^, which follows from the on-shell condition. However, to this order in 
the expansion there is a threshold of the form ~ s log s that corresponds to the non-analytic 
term that arises in supergravity in dimensional regularization around D = 10 dimensions. 

The weight-2 contribution also comes from a single Feynman diagram, in which the two 
propagators join the same pair of vertices. This gives D 2 = E 2 , so 

= 0 (7.15) 

It is notable that the vanishing of the coefficient of the 1/4-BPS interaction, is a 

special feature of D = 10 dimensions, where there are no 1/4-BPS states. The corresponding 
coefficient is non-zero for toroidal compactihcations to D = 10 — d dimensions, where 1/4- 
BPS states do exist and circulate in the loop [20,25]. 

The weight-3 contribution arises from the sum of the two contributions in (7.11). The 
expressions Di 11 = and D3 = = D3 + (/(3) (from (3.34)) lead to [4] 

S<”’‘> = |C(3) (7.16) 
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7.3 The sum of all weight-4 terms and the D^IZ'^ interaction 

The expression for Di{T) given in (5.12) enters in the sum of weight-4 terms in in 

(7.12), imply that 

jP®(r) = C2,u(r) - 1 E,(t) + 1 E^(Tf (7.17) 

Using the Laplace equation in (1.6) satished by C 2 ,i,i('r) and the one satished by the Eisen- 
stein series in (1.4), we can re-express this quantity as 

j(2-o)(r) = iAU2,i,i(r) -5E4(r) + E2{Tf (7.18) 

Now we can evaluate the integral over the fundamental domain with cutoff L, by integrating 
by part the Laplacian 

[ dfii AC2,i,i{r) - ^ = f drA^C2,i,i{T)-^E^{T)]dTi (7.19) 

\ / Jo \ / t2=L 

Using the constant term expansion for the Eisenstein series in (3.23) and for C' 2 ,i,i(r) in (5.6), 
we obtain in the large-L limit 


dfj,i ( - AU 2 ,i,i(r) —— Ai74(r) J = —C(3) -l- positive powers of L -|- 0{L ^) (7.20) 
V 2 ’ ’ 12 / 90 


The square of the Eisenstein series is integrated using Green’s theorem in the cutoff funda¬ 
mental domain leading to [4,43] 


EpE.ir>EAr) ^ <(^^0 f)- 


> Mil '^2 


+ S 5 _|_ g' _ X 

4C(2s - l)C(2s' - l)r(s - i)r(s' - i) 


(7.21) 


r(s)r(s')7r^+"'-^ 

^ 4C(25)C(25--l)r(5--^) L--' 
r(s')7r^’''^'“^ s — s' 

4a2s')a2s-l)T{s-l) L-'- ^ 
r(s) 7 r"+*'-i s-s' 


s -f- — 1 


The correction term Cl is given by the integral 


Cr = 


d/ii E ,{ t ) Es >{ t ) 


(7.22) 


r(5)r(5^ 


r 7 

-^n^+^'"Vi_25(n)o-i_2s'(n) / K^_i{27rnT2) K^,_i{27rnT2) — 

) Jr2>L " " A 
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where Fs{t) is the non-zero Fourier mode part of the Eisenstein series in the second line 
of (3.23). But bounding the Bessel function \Ks{2nx)\ < e~‘^™ / ^/x ioi x>l one gets 


0 <Cl < 


32 

r(s)r(»')i 


n>0 


^ 0 


L —)■ CXD 


(7.23) 


Specializing this formula to the case s = s' = 2 we hnd that 


J Mil 

where the scale fi is given by 


O-TT TT^ 

dhi E^irf = -C(3) log(L//r) + + 0{L-^) 


(7.24) 


, 1 , o , C'(3) C'(4) 


(7.25) 


Collecting all the contributions, the integral of j42’°^(t) over the cutoff fundamental domain 
is given by 

r TT Stt 

/ d/ii j4^’°)(r) = —C(3) + 77C(3) log(L//i) + positive powers of L + C>(L"^) (7.26) 

Jmil 90 45 

where we will drop terms that are positive powers of L (since these cancel against other terms 
arising from the low energy expansion of the non-analytic threshold terms, as described 
earlier). The logarithmic term, log(L//i) arises from the integral of the term linear in T 2 
in the large-r 2 expansion of 7^2(r)^. The logL cancels upon addition of the non-analytic 
threshold term, which was given in [25, eq. (4.41)] by, 

- log{—a'sL7r/fl4) -f (s —)■ t) -|- (s —)■ m) , (7.27) 

45 

where /i 4 is another constant. The C(3)/90 in (7.26) arises as a boundary term from the 
integral of the term involving the Laplacian. It can be absorbed into a redehnition of the 
scale, p., in log(L//i), after recalling that = (y\l2 (and remembering that the 

Mandelstam invariants have been rescaled by q;'/2 in this paper relative to [25]). 

In summary, the term in the low energy expansion of order a\ = 0(<s^) has a logarithmic 
normal threshold. An additional constant proportional to C(3) arising from the integral over 
hnite values of y was incorrectly omitted in [25]. However, the separation of the analytic 
and non-analytic contributions is ambiguous and the additive constant can be absorbed into 
the scale of the logarithm. With this dehnition, the analytic term has zero coefficient, i.e. 

2(2,0) ^ Q_ 
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7.4 The sum of all weight-5 terms and the 2 ^ 10 -^4 interaction 

In section 6 we motivated strong conjectures that express the modular functions that enter 
into in (7.13) as the sum of D 2 ,is,i{t) = C 344 (r) and sums of products of Eisenstein 

series, where satishes the inhomogeneous Laplace equation (3.19). 

To be precise, the two equations for weight 5 coefficients that correspond to Ca,h,c{T) 
functions (namely, C 2 , 2 ,i('r) = L)i^i^i^i;i(r) and = D 2 ,i,i,i{t)) have been rigorously 

derived. The other equations are based on the assumption that knowledge of all the terms 
that are powers of l/y, as well as all the power-behaved corrections to terms of order in 
the large-y limit, are sufficient to determine the full modular function. Using this information 
it follows that the total contribution of the weight-5 terms to j*'^’^^(t) is given by. 


j(i’b(^) = A ^336 6 * 3 , 1 ,i(r) -h 240^2(r) E 3 (r) -h 48U2(r) C(3) 
From (6.1) we have. 


1632 ^,, 144^,, 

-E,{t) + —C(5) 


5 


(7.28) 


C* 3 ,i,i('r) — xAUs i i(r) — —i 75 (r) 
6 15 


^-E2E,ir)-^-^ 

3 2 3V ; gQ 


(7.29) 


so that we hnd, 

^ |^56AU3,i,i(r) - ^^E5(r) -h 464E2(r) E^i^r) + 48E2(r)C(3) -h 

(7.30) 

We will now integrate this expression over the cutoff fundamental domain Ain- The integral 
of E 2 {t) E^i^r) is obtained by setting s = 2 and s' = 3 in the (7.21) resulting into 

As described earlier, terms that are powers of L cancel with corresponding terms arising 
from the low energy expansion of the non-analytic terms that come from the integral over 
AiiR. The absence of an L-independent term in (7.31) implies it does not contribute to the 
weight-5 coefficient. Likewise, the integrals of the terms proportional to E^{t) and E 2 {t) 
in (7.28) do not contribute L-independent terms. 

The integral of AC' 3 ,i,i(r) gives a boundary term proportional to dr 2 f dri 6 * 3 , 1 ,i|t- 2 =l, 
which only contains terms that behave as with a 7 ^ 0. This follows from the fact that 
there is no term proportional to T 2 in the large-r 2 expansion in (6.2). Thus, we hnd, 

[ d/ii 6 * 3 , 1 ,i(r) =-h powers of L-h C>(L"^) (7.32) 

Jmil 3 60 
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where the powers of L cancel with the integral over M-r. Substituting in (7.30) we hnd that, 


SIW'=/' d^,j<U)(r) = ^^C(5) (7.33) 

(which differs from the incorrect expression given in [25]). It will be interesting to compare 
this value with result from space-time supersymmetry and duality, when they will become 
available. 


8 Discussion and Comments 

This paper has been concerned with properties of non-holomorphic modular invariants that 
arise in the low energy expansion of the four-graviton amplitude in Type II superstring 
theory. These are dehned by Feynman diagrams for a free massless bosonic held theory on a 
torus with complex structure modulus r. Any diagram of weight w (i.e. of order («')"' in the 
low energy expansion) consists of w powers of the world-sheet Green function joining pairs 
of vertices. If the diagram has i loops it reduces to a multiple sum over the 2i components 
of the integer-valued world-sheet loop momenta. The i = 1 diagrams simply reduce to 
standard non-holomorphic Eisenstein series. We determined the exact functional form of 
the inhnite class of i = 2 invariants denoted Ca^b,c that are dehned by setting /c = 0 in (3.1) 
and have weight w = a + b + c. We have furthermore provided evidence for conjectured 
expressions for the w = 4 and w = 5 functions 714, 7 ^2, 2 , 1 , and D^, which have i = 3, 3 and 
4 loops, respectively. The conjectured expressions of these functions, which enter into the 
expansion of the four-graviton amplitude, are given in (5.12), (6.9) and (6.14). Making use 
of these expressions, we determined the integral over r of the sum of all the diagrams that 
contribute to the low energy expansion of the four-graviton amplitude at weights 4 and 5, 
thereby determining the coefficients of the one-loop ehective and interactions. 

An obvious challenge is to provide a proof of the differential and algebraic relations for the 
modular functions 7 I 4 , and 7^2, 2 , 1 - Referring to (5.12), (6.9) and (6.14) we see that any of 
these functions of weight w and ^ loops is expressed as a weight-tc polynomial of 71-functions 
with < ^ loops with rational coefficients together with powers of Riemann zeta values^®. 
These relations were motivated in sections 5 and 6 by analysis of the asymptotic behaviour 
of the various functions in the limit T 2 ^ 00 , including the hrst exponential corrections of 
order 0(e“^’^^^). This is well illustrated for the w = 5, i = 4 function that is conjectured 
to satisfy the relation (6.9). The right-hand side of this relation is the unique combination 
of 71-functions of weight w = 5 and depth 27 < 8 that matches our asymptotic analysis of 


is notable that this relationship between weight-w combinations of D functions has a structure anal¬ 
ogous to that of the relationship between multi-zeta functions, with the number of loops, £, playing the role 
of “depth”. 
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D^. Our conjecture is that F = 0, where F = — 6 OO 344 — IOF 2 O 144 + 48 F 5 — 16^(5). 

If it turns out that F 7 ^ 0, then F is a cuspidal function - a non-holomorphic modular 
function that has the behavior as r 2 —)■ 00 . This is not expressible in terms of 

any combination of the F-functions, Eisenstein series and zeta values. If, however, these 
conjectures are correct they suggest a possible extension to polynomial relations between 
F-functions of any weight. 

In addition, in order to understand the complete basis of such functions at any weight it 
is necessary to extend the considerations of this paper to consider the one-loop amplitude 
with N scattering particles, massless and/or massive. When N > 4 the Feynman rules for 
the coefficient modular functions given in section 2 are modified to include 2{N — 4) factors 
of world-sheet momenta in the numerators. Many of the resulting Feynman diagrams can 
be reduced to diagrams with fewer powers of momenta by integration by parts inside the 
Feynman diagram, as was discussed in [39] in the N = 5 and N = 6 cases. There are 
no further weight 4 diagrams beyond those that arise in the expansion of the N = 4 case 
discussed in this paper. At weight 5, however, the complete set includes two extra diagrams 
that first arise when = 5, At weight 6 the complete set includes seven new diagrams that 
Erst arise when N = 6. Understanding the relationships between the diagrams arising at 
higher weight and higher N is important for understanding the complete basis at arbitrary 
weight. 

The results of this paper generalize straightforwardly to the genus-one Type II string 
theory amplitude compactihed on a d-dimensional torus, T'^. In this case the amplitude 
(2.3) is modified by multiplying Bi in the integrand by the standard lattice factor rrf rf(r, pd), 
which ensures that Kaluza-Klein charges and string windings are acconnted for appropriately. 
The qnantities pd E SO{d, d, Z)\SO{d, d, M)/ {SO{d, M) 0 S'O(d, d, M)) parameterize the coset 
space appropriate to the T-dnality gronp SO{d, d). It is notable that the part of the genus- 
one amplitude that is non-analytic in s, t and u is independent of the compactihcation 
modnli. Consequently, the non-analytic dependence can be eliminated by consideration of 
compactihcation on tori of different modnli, pd and p'^. This provides a method of nniqnely 
specifying the local part of the genus-one effective interactions At higher genus the 

interplay between the analytic and non-analytic parts of the amplitnde is more intricate, as 
can be argned from considerations based on nnitarity. 

There are several directions in which it would be interesting to extend these results, both 
because of their mathematical interest and because of their connections to string theory. 
One direction that would obviously be of interest is to develop a complete nnderstanding 
of the structure of all the modular invariants at weight w for all valnes of w that arise in 
the expansion of the fonr-graviton amplitnde. One elegant way of establishing this might be 
by considering the action of the Laplacian directly on Fi(s,t,n|r), which is the generating 
fnnction for the modular invariants, rather than on the individnal terms in its expansion, 
generalizing the discussion in section 2 . 6 . 
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Another possible approach to determining the low energy expansion of closed-string one- 
loop amplitudes may be to exploit the relationship to the structure of the open-string low en¬ 
ergy expansion discussed in [35] in a manner analogous to the tree-level KLT relations. In the 
tree-level case the open-string expansion is a series of terms with coefficients that are rational 
multiples of multi-zeta values, which are special values of multiple poly-logarithms [26]. In 
the case of the Type II string these generalize to special values of single-valued multiple poly¬ 
logarithms that give coefficients that are (odd weight) single-valued multi-zeta values [28,30] 
The one-loop amplitude involves an elliptic generalization of these considerations. 

A different direction of interest would be to extend the analysis of this paper to the low 
energy expansion of higher-genus loop amplitudes. For genus h > 2, the world-sheets are 
really super Riemann surfaces, and the integrations are to be performed over their associated 
supermoduli spaces (for overviews see [44,45]). For genus h > 3, it may not be possible to 
project this integration onto an integration over the corresponding bosonic moduli space [46]. 
Although the modular group Sp{2h, Z) for higher genus h > 3 still exists and must play a 
leading role, it is unclear at this time which role, if any, will be played by modular forms 
under this group. 

However, for genus h = 2, higher derivative interactions involving only NS-NS fields may 
be reduced to integrals of an Sp{4, Z) modular invariant over the moduli space of compact 
genus 2 Riemann surfaces [47, 48]. In this case the leading term in the low energy limit 
is D^TZ^, which has a coefficient that is simply a constant (which was determined in [8]) 
integrated over genus-two moduli space. It is therefore proportional to the volume of the 
Siegel fundamental domain (which coincides with moduli space for h = 1,2,3). The next 
term in the low energy expansion is the D^TZ'^ interaction, which has a coefficient with 
a 5*^(4, Z,)-invariant integrand that is linear in the genus-two Green function. This was 
shown [24] to be proportional to the first power of the genus-two invariant (p introduced by 
Zhang [49] and Kawazumi [50]. The evaluation of the integral over genus-two moduli space 
in [12] was made possible by the realization that the integrand p satishes an eigenvalue 
equation Ap = 5p in the interior of moduli space, where A is the 5*^(4, Z)-invariant Laplace- 
Beltrami operator on the genus-two moduli space. Our hope is that the pattern of Laplace 
equations satished by the genus-one coefficients obtained in this paper, hints at an analogous 
pattern at at genus two and at higher genus. 
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A Asymptotic expansion of by integration 

The asymptotic expansion near the cusp as y = 71 T 2 tends to 00 of the modular functions 
which occur in this paper may be evaluated with the help of a number of different approaches. 
The hrst method relies on carrying out the integrations over the world-sheet of combinations 
of the propagator Q, and subsequent evaluation of the remaining sums. This method has 
been used already extensively in [25] to evaluate the pure power part of the asymptotic 
expansion. In this appendix, we shall extend the method to evaluating also the leading 
exponential correction terms in the asymptotic expansion of D 4 and ZI 5 . In each case, only 
a single integration over the world-sheet is required. 

For the function 712,2 ,1 two integrations over the world-sheet would be required, and 
this makes the integration method less suitable for the calculation of the asymptotics of 
112 , 2,1 • Therefore, in the next appendix, the asymptotics of 712,2 ,1 will be calculated by direct 
evaluation of the sums over the world-sheet momenta, using the decomposition into partial 
fractions of the denominators entering the sums. 

A.l Asymptotic expansion by integration 

In a hrst approach, we shall integrate combinations of the Green function Q of (2.15) over 
the world-sheet torus S with complex structure modulus r. This method is suitable for the 
calculation of the asymptotics of the functions D^, which is given by, 

De{r) = f ^dz AdzG{z\TY (A.l) 

7s 2 

and we shall describe the method specihcally for this case. We shall parametrize the torus 
S of modulus r = Ti + iT 2 by real coordinates a,/3 E M/Z with z = a \ (3t, and use a 
representation of the Green function in which Q is expanded in Fourier modes with respect 
to a. It will be convenient to make a specihc choice for the range of jd given by |/9| < 7. We 
then have the following representation for G, 

Q{z\t) = Qi{z\t) + Q2{z\t) + Qz{z\t) (A. 2 ) 

where the functions in this decomposition are given as follows, 

Gi{z\t) = 2'KT2 - \G\ + 
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Q2{,z\t) 


y Aexp 

^ \m\ 

m^O 

m^O fc^O 


a + /3 ti) — 2TTT2\ml3\ 

a + /3 ti + kri) — 2 TiT 2 \m{k + /9)| 


(A.3) 


The above decomposition holds for |/9| < | and may be extended beyond this domain by 
periodicity in {3 with period 1. For the evaluation of such extension will not be needed, 
but it will enter for other types of diagrams. The integral over S of Qi vanishes. The function 
^3 is bounded uniformly in by an exponential in y, 


\Qz{z\T)\<C{y)e-y (A.4) 

where C{y) is power-behaved and independent of ; 2 . Therefore, will not contribute to the 
pure power part of the asymptotic expansion, while the leading exponential correction 

and its complex conjugate will involve Q-^ to first and second order. 


A.2 Evaluating the pure power part 

We begin by reviewing the calculation of the pure power part of the expansion, following [25]. 
The pure power terms are obtained by retaining the contributions from Qi and Q 2 only, and 
extracting the pure power part of the corresponding integral, 

= j da J ^ d/3 (Gi{z\t) + Q 2 {z\r)'J + C>(e"^) (A.5) 

while omitting any exponential terms that still arise in this expression. It will be useful to 
expand the integrand as follows, 

E (A.6) 

^1/2=0 

where is the Kronecker <5, while the expansion coefficients are, 

= I da y*' dl3g,{z\TY^ g2{z\rY^ + 0{e-^) (A.7) 

Carrying out hrst the integration over a gives d^^’^\ii, 1;?/) = 0 for all ii, and also readily 
produces the contribution with £2 = 0 , 

0; v) = 2^b/?6,(^|T))'- = (|)'‘ 2 U ( 1 . -<i; |) (A.8) 


59 



The hypergeometric function 2 T 1 evaluates to a rational number. For £2 > 2, we have, 


da Q 2 {z\tY^ = 




m,0 


^ |mi... 


,-2yMm 


where we have used the following notation, 


(A.9) 


m = mi + m 2 + ... + 771^2 

M = |mi| + |m2| + ... + |m£2| (A.10) 


In the remaining integration over /9, the integrand is even in fd, so we shall restriction the 
range to 0 < /9 < | and include a factor of 2. Since M > 2, the integral over (3 from | to 
cx) is uniformly bounded by a multiple of so that the pure power part of the expansion 
precisely corresponds to taking the integration domain to be 0 < /9 < cxd. The corresponding 
integral over (3 is then given by a polynomial in of degree 2^l + 1, whose coefficients 

will be denoted by P(f'i,n) and are given as follows. 


6 


Collecting all contributions, we hnd. 




n=l 




(A.ll) 


2£i+l 


S^’^\ii,i2;y) = 2j2 P{h,n)S{i 2 ,n) {2y) 


ii—n 


n=l 


where the mnltiple sum S is dehned by 


S(( 2 .n)= 




m.O 




\mi... miY 


(A.12) 


(A.13) 


These multiple sums may be expressed in terms of multi-zeta functions (see [25] and [38]), 
which are functions of complex variables si,... ,Sp defined by. 


C(Si,S2,...,Sp) = ^ 


ni>...>np>l 


Ui . . . Un 


(A.14) 


and we have the following relation, 

‘S(72,n) = 7|- 5^ 2^'’+‘^“'"C(n + 2,o,,...,a,) 


(A.15) 


ai,...,apG{l,2} 
a-j^ +... + a p = 2 ~ 2 


For special valnes of p and of the argument, the multi-zeta functions may be expressed in 
terms of polynomials in the ordinary (C-fnnctions. 
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A.3 Evaluating the leading exponential correction 

The contribution to Di{t) arising from terms in Qi and Q 2 only (while omitting ^3) do 
produce exponential corrections which are independent of ri, as may be seen explicitly 
from the fact that the a-integration over is independent of Ti. As a result, those parts 
contribute exponential corrections vf’^\y) only for k = k, and will not contribute to the 
leading exponential corrections which have k + k = l with ki, k 2 > 0. Therefore, the leading 
exponential corrections will arise from expanding in powers of Q 3 , up to order 2 included, in 
view of the bound (A.4), and we have. 


- X] <^4+£2/ 

^1/2=0 





(A.I6) 


where d^^’^\ii,i 2 ;y) is the contribution proportional to q (indicated by \q below) in. 


d^^’^\ii,i 2 ;y) = da j ^dpOiizlrY^ {Q 2 {z\t) + G-i{z\r)^ 
By inspection of the integration over a, it is clear that we have. 


«2 


(A.17) 


(A.18) 


for all £1. For all £2 > 3, the term of order is multiplied by at least one power of G2, 
and this suppresses the exponential by at least the order so that such terms do not 

contribute to the leading exponential corrections. Therefore, the contribution for £2 ^ 3 
takes the form. 


d^^'^\h,h]y) = h j j ^df3gi{z\TY^ g2{z\rY^ ^Gsizlr) 

while for £2 = 2 , we have, 

S^’^\ii,2-,y) = 2/" da [ d/3 GYzItY"- g 2 {z\T) g 3 {z\T) 


1 

' 2 
1 
2 


+ I da I d/S ggzlrY^ Gsizlr) 


(A.19) 


(A.20) 


The evaluation of these integrals proceeds along the same lines as for the pure power part, 
and makes use again of the formula (A. 11), and we have for £2 ^ 3, 


2^1+1 

vY'^\ii,i 2 ;y) = P(^i,n) S^^\l 2 ,n) {2yY^-- (A.21) 

n=\ 
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where the sums are gives as follows for £2 > 3, 




™.(^^-2)' 

with m and M dehned by (A. 10). Finally, for the special case £2 = 2, we have, 




2^1 + 1 


£1—n 


(A.22) 


(A.23) 


n=l 


Note that the sums S^^\i 2 ,n) diverge for £2 = 2 because the combination (M — 2)"' is then 
forced to vanish. In the correct formula for the case £2 = 2 given above, the part of S^^\£ 2 , n) 
that involves only the combination M” survives and gives rise to the factor ^ in (A.23). 

We shall need to evaluate these sums up to £2 = 5 included to calculate the expansions 
up to ZI5. It is straightforward to establish the following recursion relations, 

SW(3,n) = is<‘l(3,n-2) + |;(l-C(n)) 

1 12 / \ 
^«(4,n) = -A«(4,n-2) + - UC(n,l)-C(n + 2)-2n-l + 2^CWj (A.24) 

along with the initial conditions, 

^(b(3^0) = 8 ^4)(4^0) = 24 + 47r2 

^4)(3,i) = 2 = 12_7r2 + 6C(3) (A.25) 

For £2 = 5, we shall only need a single sum, given by, 

^4)(5,i) = 72 +—-48C(3) (A.26) 

-LO 

From these data, the leading exponential corrections to ZI4 and may now be readily 
obtained, and one hnds the results given respectively in (5.9) and (6.6). 


B Asymptotic expansion of i^2,2,i by summation 


One key ingredient in /l2,2,i is the appearance of a subgraph which quantum held theorists 
refer to as the one-loop vacuum polarization diagram, and which in terms of discrete world- 
sheet momenta is given by the following sum over internal loop momenta (m,n). 


T{M,N\t) 


(m,n)7^(0,0), 


TT^Imr -|- np |(m -|- M)t + n + Wp 


(B.l) 
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From its definition, we have T{—M, —N\t) = T{M, N\t) and T{M, N\t) > 0 for all M,N. 
The summation over the integer n may be carried out explicitly by decomposing the sum¬ 
mand in partial fractions with respect to the variable n, and using the following elementary 
summation formula, 


E 


1 

mr + n 


= —in 


l + q^ 

l-qm 


q = e 


2'KiT 


(B.2) 


This sum is conditionally convergent, and may be naturally dehned by taking the limit as 
A —)■ cx) of the hnite sum obtained by restricting the range of n to —A < n < A, following 
Eisenstein. When M = 0 the sum reduces to. 


r(o,iv) = 






n'^Ti^iN + nY 

nj^0,—N ' '' 


EE 

m^O n 


To 


vr^lmr -|- nP |mr + n + N\‘^ 


(B.3) 


In particular, we have T(0, 0) = E 2 . When 7^ 0, we decompose into partial fractions in n 
and carry out the sums over n, and we hnd. 


r(o,iv) = 


2r| 


6r| 


E 


2to 


1 + g" 


c.c. 


?)N‘^ n'^N'^ ~ 9 

When M 7^ 0 the variable N is unrestricted, and the summation over n gives, 

iT2 iT2 


(B.4) 


T{M,N) = 


2r| 

3|Mr + A^P \n^M{MT + N)^ n^M{Mf + Nf 
T 2 1 1 + q 


E 

m^O 


-I- C.C. 


nm {Mr + N){Mt -|- A^ — 2imT2) I — q'‘ 

Note that both series in m are absolutely convergent for y > 0. We have for example. 


(B.5) 


Do = 


E 


T2 


D 


(M,Af)^(0,0) 

T{M,N) 

{M,N) 

2,2,1 = E 


7r|Mr -|- A^p 
2 


r(M, N) 


T2 


(M,V)^(0,0) 


7r|Mr -|- A^|- 


r(M, Nf 


(B.6) 


In terms of T(M, A^), the summation needed for will involve two internal momenta, which 
renders this approach less effective for calculating the asymptotic expansion of D^. 
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B.l The power and leading exponential approximations 

As we shall be mostly interested for the moment in computing the pure power and leading 
order exponential corrections, we may right away approximate the two-point function to this 
order. Subsequent summations over the “external momenta” M, N will not interfere with 
the approximation. We write the following asymptotic expansion, 

r(M, N) = To{M, N) + q Ti{M, N) + qT^{M,N) + 0{\q\^) (B.7) 


Note that Ti{M,N) need not be real. These leading components may be read off from the 
formulas obtained above. For the special cases where M = 0, we hnd. 


ro(o,o) 

71(0,0) 


45 y 
1 

2 + - 
y 


as well as when N ^ 0, 


ro(o,iv) 

ri(o,iv) 


2 r| 

3iV2 


6 r| 4 r 2 1 

+ 4m^r|) 


At2 

+ 4r|) 


For the case M 7 ^ 0, we obtain. 


To{M,N) 


r,{M,N) 


2 r| 


iT2 


IT2 


3|Mr + iVp + N)^ Ti‘^M{Mf + Nf 

T2 


E 

m^O 

2r2 


[7r|m|(Mr -|- N){Mt -|- TV — 2imT2) 

1 1 

+ 


-I- c.c. 


7 r(Mr -f N) \Mt + N — 2iT2 Mr -|- TV -|- 2iT2 


(B.8) 


(B.9) 


(B.IO) 


B.1.1 Useful summation formulas 

When calculating the sums needed to evaluate the C- and U-functions, it is very convenient 
to organize the calculation in terms of the poly-logarithms dehned by, 

'^{n,z) = —lnr{z) (B.ll) 
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where ^(;2) is often used for \l/(0,;2). Finite sums are conveniently expressed in terms of 
these functions, since we have, 

m—1 ^ 

'F(n, m) - 'F(n, 1) = (-)'^r(n + 1) ^ 

i=l 

with the following explicit evaluations of their values at z = 1, 

'F(0,1) = -7 

^(n,l) = (-)"+ir(n + l)C(n+l) 

The following sums may be recast in terms of multi-zeta functions, 

m=l 

where convergence requires k > 2. A further generalization of this result will be needed for 
sums quadratic in T, and we have, 

OO 

m=l a=l,2 

= (-)^i+^2r(A;i)r(A;2) (({k, ki + k2) + C(A:, k2) + ({k, k^, fci)) (B.15) 

Finally, the following two sums will also enter, 

= ^C(7)-16C(2)C(5)-4C(3)C(4) 

m=l 

^ 4.(1,2,n.)-4.(1,!) _ _ 4<;(3)<;(4) (B.16) 

mP 1 

m=l 

Although evaluations of these sums have been claimed in [51], the expressions given for them 
in formulas (1.10) and (1-13) of that paper do not stand the test of numerical verihcation. 
The validity of the above formulas (B.16) has then been verified numerically up to a relative 
precision of 10“^® using MAPLE. 

B.2 Evaluating the pure power part ^^ 2 , 2,1 

We decompose the calculation of £*2,2,1 into the contributions coming from M = 0 and those 
coming from M 7^ 0, 

D2,2,i = B'°’ + B*'* (B.17) 


(B.12) 

(B43) 

(B.14) 
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where each piece is given by, 


®'“’ = T.i^.no,N) 




TfiV^ 


®'‘’ = EE 


T2 


M^O N 

We shall now evaluate each contribution in turn. 


tt\Mt + iVp 


r(M, Nf 


(B.18) 


B.2.1 Evaluating the power terms in 2)^^^ 

Retaining only the power terms, the contribution itself splits into three parts, according 
to whether it involves one, two, or three sums, 

S)(o) = + 2 )® (B.19) 


where these terms are the power parts of the following sums. 


^(0) ^ 'sp ^2 ( _ 6r| \ ^ 

1 ^vriV^VdiV^ TT^iVV 

^( 0 ) ^ gV ^2 /2r| 6r| \ 2 r 2 1 

^ ySiV^ TT^iV^/ Tim (iV^ + 4m^r|) 

CXD CXD \0 i 

2>(0) ^ ^2 ’sp -sp (2r2) _4_ 

^ ^ + 4/c2r|)(iV2 + 4m^r|) 


The summations over N are carried out with the help of (B.2). The remaining summations 
over k and m are then transformed into sums that give rise to ^-functions and their products. 
The results are as follows. 


^(0) ^ 8y^ 

^ 467775 

2)(°) = + 

2 945 ^^^ 45 5670 |/ 3y^ 25200 16y^ 

®7 = ^ + 4(2C(7)-<;(3)C(4)-C(2)C(5)) + diph (B,21) 


B.2.2 Evaluating the power terms in 

Retaining only the power terms, the contribution 2 )^) itself splits into three parts, according 
to whether it involves one, two, or three sums in addition to the sum over M 7 ^ 0 , 

2 )(i) = 2)^^^ + (B.22) 
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where 


2)« 


M^O N 


T2 


2rn 


7r|Mr + 7V|2 \3\Mt + N\^ { + N)^ 


IT2 


2 EE 


T2 


M^O N 

xE 

m^O 


2r^ 


*7-2 


7r|Mr + iV|2 \3\Mt + N\^ { tt^M{Mt + N)^ 

T2 


+ c.c. 


+ c.c. 


7r|m|(Mr + N){Mt + N — 2imT2] 


+ c.c. 


EE 


T2 


E 


T2 


M^o N t^Wt + N\^ ^^\Ti\k\{MT + N){Mt + N - 2ikT2) 


c.c. 


xE 

m^O 


T2 


7r|m|(Mr + N){Mt + N — 2imT2] 


+ c.c. 


(B.23) 


In each case, the sum over N is computed hrst with the help of (B.2). For 2)^^^ the power 
contributions to the remaining sum over M 7 ^ 0 are computed in terms of ^-functions. For 
5 ) 2 ^^ one needs to appeal to a few simple multi-zeta functions. For the sum over N 
reveals that its entire contribution is proportional to y~^. 

In summary, the results are as follows. 


(1) ^ C(5) 5C(7) 21C(9) 
^ 3 6|/2 lQy4: 


where is a number which is independent of r. Collecting all the contributions to the value 
of F* 2 , 2 ,i up to the knowledge of B, which affects only a single power of y, we hnd. 


V 


( 0 , 0 ) 

2,2,1 


(y) 


8y^ V .. 3 ^ 13C(5) _ C(3)C(5) ^ 

467775 945 ^^ ^ 45 y^ 8y^ 

+ (^ + 40C(7) - 16C(3)C(4) - 16C(2)C(5)) 


(B.25) 


The contribution B is the most challenging to compute. It is given by the power part of the 
following sum. 


B 


16'^! 1 f f _ ]_ _ 

^ ^,^\MT + N\^\{^^\\k\{MT + N){Mf + N 



(B.26) 
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Within the power part approximation, is a constant, as becomes evident upon carrying 
out the summation over N using (B.2). To carry out the summation over N, we decompose 
into partial fractions. The order of the poles jumps when M = —k, and these contributions 
need to be isolated, and handled separately. We shall omit the details of the calculation 
here, and simply quote the result, 

B = -26C(7) + 16C(3)C(4) + 16C(2)C(5) (B.27) 

which yields the result for T>^ 2 i{y) quoted in ( 6 . 12 ). 


B.3 Evaluating the leading exponential terms 

The leading exponential contribution to splits into three parts, according to whether it 
involves zero, one or two sums, and we shall denote these contributions by S, 


where 








= 0 

^ ^ ~ ^ 7r(iV2 + V) 

c(o) = ^ _ 

^ ^ + Ay'^){N‘^ + Arn^y"^) 


(B.28) 


(B.29) 


The summations over N are carried out with the help of the formulas of appendix A. The 
remaining summation over m are then transformed into sums that give rise to ^-functions. 
The results are as follows. 




’( 0 ) 

1 

= 0 

’( 0 ) 

Ay^ 

'2 

945 

’( 0 ) 

_ A. 

3 

3|/^ 


er = 


2 1 
45 &y 


1 3 




y" 


TT 

6|/^ 


TT 


+ ^C(3) + —C(5) 


90|/2 ^ 2|/3 


(B.30) 


B.3.1 Evaluating the leading exponential terms in pB) 

The leading exponential contribution to splits into three parts, according to whether it 
involves zero, one, or two, sums in addition to the sum over M 7 ^ 0 , 


p{i) 


= q[£l^Ys? +8^^^ 
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9 


(B.31) 



where 



M^O N 


y 

7 r|Mr + iVp 


ro{M,N)ri{M,N) 


(B.32) 


The remaining sums are the parts linear in q of the sums dehned earlier, 

2y^ 


= EE 


y 


N 


7i\Mt + iV|2 \3\Mt + iV|2 [ + iV)3 


w 


9d” 


2 EE 


y 


MjtO N 


2y^ 


7i\Mt + iV|2 \3\Mt + iV|2 [ + iV)3 


w 


“h C.C. 


+ C.C. 


xE 

m^O 


y 


7 r|m|(Mr + N){Mt + N — 2imy) 


+ C.C. 


^4" = EE 


E 


M^o N + N\^ ^^\Ti\k\{MT + N){Mt + N - 2iky) 


+ C.C. 


xE 

mT^O 


y 


7 r|m|(Mr + N){Mt + N — 2imy) 


C.C. 


(B.33) 


In each case, the sums over N are computed hrst, with the help of the formulas of appendix 
A. The contributions involving the sums over k and m are further partitioned according to 
whether either A; or m or both are equal to —M, or not. 


y ( 1 ) 

_ y(l) 1 y(l) 1 y(l) 

— ^Oa “T ^Ob ‘xQc 


y (1) 

y(l) _i_ y(l) 

~ ^2a “T ^2fe 


y (1) 

y(l) _i_ y(l) _i_ y(l) 

~ ^3a “T ^3b “T ‘^3c 

(B.34) 


Here is the contribution from the part of %{M, N) that does not involve a sum over m, 
while 8qI^ is the contribution that arises from the sum part of %{M,N) of the single term 
m = —M, and hnally 8^]} is the contribution from the sum part in %{M, N) with M ^ —m. 
Furthermore, 8^^^ is the contribution from m = —M while 8^1^ is the contribution from 
m 7 ^ —M. Similarly, £^ 3 ^^ is the contribution from k = m = —M, while £^ 3 ^^ is the combined 
contribution from k = —M, m 7 ^ —M and k 7 ^ —M, m = —M, while 8^1^ is the contribution 
from k,m ^ —M. These contributions are found to be given by, 

1 3 C(5) 

3y 4|/3 y3 
17 TT^ TT^ 

V ~ ^ ^ 9(V 


y ( 1 ) _ 

^Oa ~ 


y (1) _ 

^Ob ~ 
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c(l) _ 

^Oc ~ 


= 


c(l) 

^2a 


c(l) _ 

^2b ~ 

c(l) _ 

‘"3a ~ 


c(l) _ 

‘"36 ~ 

c(l) _ 

^3c — 


21 TT^ C(3) 
+ H-^ 


%2 2|/2 




68y^ , 122 |/ ^ 8 
^ ^ ^ 3 ^ 


8 ^ 37 ^ 21 ^ 21 

% ^ IV ^ ^ IV 


128p^ 16 19 11 11 

45 ^ 3 3y 2|/2 4^3 


^ 22 14C(3) 6 

3y 3y y"^ 

4y^ Ay ^ 7 7 


4C(3) 11 

o ' 


3C(3) C(5) 


V 


Ay^ 2y^ 


2y^ 


_4-6+!CM 

y y 

1 

V 


2y 


7 , 2C(3) 

~r o 




(B.35) 


The sum of these contributions yields (6.15). 


C Preliminary numerical study of modular functions 

The analysis of the properties of the Ca^^c functions in section 3 and the asymptotic properties 
of certain D and C functions in the y ^ oo limit in sections 5 and 6 has led us to a number of 
conjectured relations between these functions. Although we have a complete understanding 
of the relationships between the Ca,b,c functions at any weight w, our analysis of other weight- 
w functions has been limited to conjectured relationships for functions with w < 5. We do 
not have a rigorous mathematical procedure for establishing the validity of these conjectures, 
beyond the striking agreement of many terms in the asymptotic expansions. In the absence 
of an analytic method we have approached the problem of establishing such relationships by 
numerical methods. 

Even though at present our attempts are very primitive we believe they demonstrate the 
possibility of using numerical approximations to the modular functions to sufficient accuracy 
to test the validity of our conjectures. The conjectured relations that we want to test are 
summarized by: 

• One relation at weight 4, which we recall here from (1.10), 

C*i,1,1,1 = 24(7*2,1,1 + 371/2 + I 8 E 4 (C/.l) 

We also recall the alternative notation, (7*1,1,1,1 = D 4 and (7*2,1,1 = 77*2,1,i- 

• Three relations at weight 5, which we recall here from (1.11), 

406*2,1,1,1 = 3006*3,1,1 + I2OE2E3 - 276^5 + 7C(5) 
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(C.2) 


~ 60C*3^i^i + 10-E/2-E'3 ~ 48-£/5 + 10(^(3)-E/2 + 16C(5) 
101^2^2,1 = 20(73^1^1 — 4ii^5 + 3C(5) 


In the alternative notation we have (^2,1,i,i = -^3,1,1, C's,!,! = -D2,1,1,1, Ci,1,1,1,1 = D^. 

Our numerical procedure is based on evaluating the expressions for the D functions by 
using the integral representation in terms of the Green function, 


-^^12,•••/34 (^) 



(C.3) 


where G{z\t) was defined in (2.15). 

We have obtained the best compromise between precision and running time by using the 
Divonne implementation for Mathematica described in [52].^® 

The asymptotic expansions for large T 2 described earlier include power behaved correc¬ 
tions to the first exponentially suppressed term proportional to and ignore corrections 

of order which are highly suppressed for relatively small values of T 2 . Therefore, in 

order to provide a numerical test that goes beyond the asymptotic expansions which mo¬ 
tivated the conjectured relationships, we need to evaluate the functions with considerable 
accuracy, which we have found difficult to achieve. 

In order to keep the computation time to at most the order of a few days, in the following 
we have attempted a precision of e = 10“^, 10“^, 10“^ or 10“^ with about 500 x 10® sampling 
points, and using four values of the complex structure parameter r G {i, ^ + \ + 7i}. 


Function 

e 

T = i 

r = iy/2 

r = i + 7i 

r=i + ^^ 

Ds 

0 

1 

1.0003 

0.9995 

0.9999 

1.1851 


10-® 

0.996 

0.998 

1.009 

1.024 

D, 

10-2 

0.95 

0.97 

1.12 

0.93 

- 02 , 2,1 

10-2 

0.96 

0.92 

1.07 

0.84 

- 03 , 1,1 

10-2 

1.10 

0.87 

1.05 

1.37 


Table 1: The ratio of left-hand and right-hand side of the identity 


The results that we will now describe are summarized in Table 1 and Table 2. As a 
benchmark, we first numerically tested the weight-3 relation, 

I^3(r) = Gi,i,i(r)=E3(r) + C(3) (C.4) 

^®We would like to thank the Niels Bohr Institute and the IHES for allowing us to use their computer 
resources. 
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Function 


T = i 

r = iy/2 

r = i + 7f 


Ds 

numerical 

asymptotic 

1.35189 

1.35227 

1.42809 

1.42744 

23.7143 

23.7119 

1.33559 

1.33542 

D4 

numerical 

asymptotic 

5.59442 

5.60153 

6.118 

6.11683 

265.665 

265.577 

5.47065 

5.46992 

D5 

numerical 

asymptotic 

27.5931 

27.6018 

29.8883 

29.9017 

1330.99 

1331.95 

27.0271 

27.0225 

F*2,2,1 

numerical 

asymptotic 

0.417141 

0.419068 

0.475233 

0.463207 

90.6985 

90.7233 

0.398955 

0.409520 

-^ 3 , 1,1 

numerical 

asymptotic* 

0.701572 

0.769175 

1.05760 

1.02236 

474.804 

474.913 

0.733569 

0.778507 


Table 2: Comparison of numerical and asymptotic values 


which is known to be correct since we proved it analytically in sub-section (3.3). The results 
are shown in the first row of each table. In this case we were able to achieve e = 10“*^, as 
noted in the second column of Table 1. Columns 3-6 of Table 1 show the ratio of the left 
and right-hand sides of equation (C.4) obtained by numerical evaluation of the functions 
D^, with the help of the integral representation in (C.3). This ratio is equal to 1 to the 
stated accuracy. In the first row of Table 2 we list the numerical estimate of the value of the 
function D^{t) at each value of r, while the second row gives the value of the function as 
determined from the asymptotic estimate, which keeps the pure power terms and the first 
exponential correction in the exact expression in (C.4). The results for the four conjectured 
relations in (C.l) and (C.2) are presented in the rows of Tables 1 and 2 underneath the 
results. 


Although all the results are consistent with the conjectured relations, they are not yet 
more accurate than the asymptotic expansions that include the first exponentially suppressed 
terms that we studied analytically in the body of the text. For this reason we are not able to 
state with conviction that the numerical data conhrms the conjectured relationships beyond 
the asymptotic analysis, but we have included this appendix since it does suggest that 
numerical methods are well within reach given a little more sophistication. 


The * on the last row in Table 2 indicates the fact that in this row only the pure power 
part of the asymptotic expansion has been included, since the leading exponential correction 
has not, in fact, been calculated yet. 
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